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ABSTRACT 

The most accurate measurements of magnetic fields in star-forming gas are based on 
the Zeeman observations analyzed by Crutcher et al. (2010). We show that their finding 
that the 3D magnetic field scales approximately as density** ®^ can also be obtained 
from analysis of the observed line-of-sight fields. We present two large-scale AMR 
MHD simulations of several thousand Mq of turbulent, isothermal, self-gravitating 
gas, one with a strong initial magnetic field (Alfven Mach number A4a,o = 1) and 
one with a weak initial field (AIa.o = 10). We construct samples of the 100 most 
massive clumps in each simulation and show that they exhibit a power-law relation 
between field strength and density (hn) in excellent agreement with the observed one. 
Our results imply that the average field in molecular clumps in the interstellar medium 
is (i3tot(hH)) ~ 42hH®4 pG. Furthermore, the median value of the ratio of the line-of- 
sight field to density®'®® in the simulations is within a factor of about (1.3, 1.7) of the 
observed value for the strong and weak field cases, respectively. The median value of 
the mass-to-flux ratio, normalized to the critical value, is 70 per cent of the line-of-sight 
value. This is larger than the 50 per cent usually cited for spherical clouds because the 
actual mass-to-flux ratio depends on the volume-weighted field, whereas the observed 
one depends on the mass-weighted field. Our results indicate that the typical molecular 
clump in the ISM is significantly supercritical (~ factor of 3). The results of our 
strong-field model are in very good quantitative agreement with the observations of 
Li et al. (2009), which show a strong correlation in field orientation between small 
and large scales. Because there is a negligible correlation in the weak-field model, we 
conclude that molecular clouds form from strongly magnetized (although magnetically 
supercritical) gas, in agreement with the conclusion of Li et al. (2009). 

Key -words: Magnetic fields—MHD—ISM: magnetic fields—ISM: kinematics and 
dynamics—stars:formation 


1 INTRODUCTION 

Since the detection of magnetic fields in the interstellar 
medium (ISM), they have played a central role in star for¬ 
mation theory. The first estimate of the strength of the in¬ 
terstellar field was made by Leverett Davis (Davis 1951), 
who inferred a field strength of 10“"^ G based on a compari¬ 
son of the observed fluctuations in the orientation of optical 
polarization vectors and the observed amplitude of turbu¬ 
lent motions in the ISM. Two years later, Chandrasekhar 
& Fermi (1953) used the same method (with the correc¬ 
tion that it is the ID velocity dispersion rather than the 
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3D value that enters) with different values for the density 
and velocity dispersion to obtain B = 7.2 pG, remarkably 
close to the current best estimate for the field strength in 
the diffuse ISM B = 6.0 ± 1.8 pG (Heiles & Troland 2005). 
(More precisely, Chandrasekhar & Fermi determined the 
field in the plane of the sky, so their value for the total field 
strength would be 8.3 pG.) This method of determining the 
field strength is generally referred to as the Chandrasekhar- 
Fermi method, but in view of the fact that it was originally 
and independently discovered by Davis, as acknowledged by 
Chandrasekhar & Fermi (1953), we shall refer to it as the 
Davis-Chandrasekhar-Fermi (DCF) method. 

In a classic paper, Mestel & Spitzer (1956) analyzed 
the implications of interstellar magnetic fields for star for¬ 
mation. They pointed out that the magnetic energy scales 
with the gravitational energy, so that the magnetic field sets 
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a lower limit on the mass of a gas cloud that can undergo 
gravitational collapse. In modern terminology, they discov¬ 
ered that there is a magnetic critical mass; masses above this 
can undergo gravitational collapse (the magnetically super¬ 
critical case), whereas lower masses cannot (the magneti¬ 
cally subcritical case). They estimated that the magnetic 
critical mass was at least 500 M©. Subsequent observations 
have shown that diffuse HI clouds (the only type of inter¬ 
stellar cloud known to Mestel and Spitzer) are indeed mag¬ 
netically sub-critical (Heiles 1997; Heiles & Troland 2004). 
Since 500 M© is much larger than the mass of a star, Mestel 
& Spitzer (1956) argued that star formation occurred via 
ambipolar diffusion, in which the neutrals contracted rela¬ 
tive to the ions that were tied to the magnetic field. They 
further speculated that the field might be strong enough 
to prevent the formation of a disc around the protostar, 
thereby anticipating what is informally known as the ‘pro- 
tostellar accretion-disc crisis’ decades before such discs were 
discovered observationally! 

Whereas Mestel & Spitzer (1956) argued that ambipo¬ 
lar diffusion was rapid, subsequent workers have found that 
it is about an order of magnitude slower than free-fall col¬ 
lapse (Mouschovias 1987; Shu et al. 1987). The classical pic¬ 
ture of low-mass star formation was based on the assump¬ 
tion that the molecular cloud cores out of which low-mass 
stars form are magnetically subcritical and that gravita¬ 
tional collapse is enabled by quasi-static ambipolar diffusion. 
This model has been challenged by Zeeman observations of 
dense molecular cloud clumps inside giant molecular clouds 
(GMCs), which show that on average the clumps are mag¬ 
netically supercritical-!.e., that gravity dominates-by a fac¬ 
tor of order 2 (Crutcher 1999; Falgarone et al. 2008; Troland 
& Crutcher 2008; Crutcher et al. 2010). Zeeman observations 
measure the net field along the line of sight, and that field 
is almost never strong enough to render the clumps magnet¬ 
ically subcritical (Crutcher 2012). The only clear counter 
example that we are aware of comes from DCF observa¬ 
tions of a quiescent molecular cloud, which indicated that 
the clumps in that cloud were uniformly slightly subcriti¬ 
cal (Marchwinski et al. 2012). The significant uncertainties 
associated with the DCF method are reviewed by Crutcher 
( 2012 ). 

Although gravitationally bound CMCs should be mag¬ 
netically supercritical (McKee 1989), fragmentation along 
the field lines should result in magnetically subcritical 
clumps; indeed, it is for this reason that Mestel & Spitzer 
(1956) invoked ambipolar diffusion to enable star forma¬ 
tion. The most likely reason that the clumps inside CMCs 
are not magnetically subcritical is due to the effects of the 
supersonic turbulence observed inside CMCs. Turbulence 
can accelerate ambipolar diffusion (Fatuzzo & Adams 2002; 
Zweibel 2002; Li et al. 2012b). Furthermore, Lazarian and 
his collaborators (e.g. Lazarian et al. 2012; ?, 2015) have sug¬ 
gested that magnetic reconnection in a turbulent medium (a 
process they term ‘reconnection diffusion’) leads to rapid dif¬ 
fusion of the magnetic field independent of ambipolar drift; 
this process is most effective on length scales at which the 
turbulence is super-Alfvenic . 

Large scale numerical simulations provide a powerful 
tool for interpreting the observations of magnetic fields in 
molecular clouds. Whereas observations of fields are limited 
to the component along the line of sight for Zeeman obser¬ 


vations and to the components in the plane of the sky for 
DCF observations, the full 3D field is accessible to numerical 
simulations. The majority of simulations used to study mag¬ 
netic fields in molecular clouds are ideal MHD (although see 
Oishi & Mac Low 2006; Li et al. 2008; McKee et al. 2010; 
Vazquez-Semadeni et al. 2011; Li et al. 2012b) and therefore 
do not include the acceleration of ambipolar diffusion by tur¬ 
bulence. On the other hand, such simulations can simulate 
reconnection in a turbulent medium insofar as numerical re¬ 
connection mimics the effects of actual reconnection. 

A number of simulations have been carried out to ex¬ 
plore super-Alfvenic models of molecular clouds, in which 
the energy in turbulence greatly exceeds the energy in the 
initial, uniform magnetic field (Padoan & Nordlund 1999; 
Padoan et al. 2004; Lunttila et al. 2008, 2009; Collins et al. 
2011). These simulations did not include self-gravity except 
that of Collins et al. (2011), who considered a cloud that 
was initially very magnetically supercritical. Comparison of 
the distribution of the line-of-sight magnetic field. Bios, ver¬ 
sus the column density of cloud clumps with observation 
showed good agreement. Both Lunttila et al. (2009) and 
Collins et al. (2011) use models with very large initial val¬ 
ues of the plasma beta, p = SivpCs/K = 22.2, corresponding 
to very weak magnetic helds. The turbulence amplified the 
field, however, eventually reducing /?(Brms) to about 0.2 in 
the simulation by Collins et al. (2011). The key observa¬ 
tional distinction between weak and strong mean fields is 
in the field structure: The field in the weak mean-field case 
is very tangled, whereas that in the strong mean-held case 
maintains a correlation in the orientation on large and small 
scales (Ostriker et al. 2001). Observations by Li et al. (2009) 
of the the orientations of the held on large and small scales 
led them to conclude that the mean held is strong-even an 
Alfvenic Mach number A4a = 2 is too large to be consis¬ 
tent with the data. One of the goals of our work is to revisit 
the issue of the orientation correlation in simulations that 
include self gravity. 

In this paper, we describe the results of two large- 
scale adaptive mesh rehnement (AMR) simulations with the 
ORION2 code (Li et al. 2012a) of supersonically turbulent, 
magnetized, self-gravitating molecular clouds (Li et al. 2015, 
in preparation). There are two objectives of our study: First, 
we wish to compare the results of our simulations with the 
inferred properties of magnetic helds in molecular clouds 
by Crutcher et al. (2010). Second, we wish to examine the 
magnetic properties of molecular-cloud clumps under strong 
and weak mean-held initial conditions in order to determine 
which agrees better with observation. In Section 2, we re¬ 
view the results of the Bayesian statistical analysis of the 
Zeeman observations of magnetic held strengths by Crutcher 
et al. (2010). Section 3 presents new results from an anal¬ 
ysis of this data. We describe the simulations, the scaling 
to physical variables, and the model parameters in Section 
4, and the selection of the samples of simulated clumps for 
comparison with observations in Section 5. We compare the 
line-of-sight properties of clumps from our simulations with 
those from the 68 OH and CN sources cataloged in Crutcher 
et al. (2010) in Section 6. In Section 7, we study the relation 
between the observed line-of-sight value of the mass-to-hux 
ratio with the intrinsic 3D value and describe its evolution. 
The long-standing question of whether the mean magnetic 
fields in molecular clouds is strong or weak is discussed in 


© RAS, MNRAS 000, 1-30 



Magnetized Interstellar Molecular Clouds. I. 3 


Table 1. Glossary 


Parameter^ 

Definition 

Reference 

Bo 

Magnitude of initial magnetic field: 31.6 (strong) 3.16/i.G (weak) 

Sec. 4.2 

^los 

Magnitude of mass-weighted line-of-sight field in central Gaussian beam 

Eq. (1) 

l-^lmid, z' 

Average magnitude of field in central of clump midplane 

Sec. 7 

r> 

-^mid, z' 

Magnitude of average axial field in central clump midplane 

Sec. 7 

Btot 

Total mass-weighted field, (B? H- B? -h B? )^/^ 

® ' los, a: los, y los, 

Eq. (22) 

Bv 

Magnitude of the field averaged over clump volume, {\B\)y 

Sec. 8.1 

-^rms, V 

RMS field in clump volume, ((B^}^)^/^ 

Sec. 6.2 

Cs 

Mass-weighted isothermal sound speed 

Eq. (A5) 

CA 

= 1 (box), 7v/4 (sphere) 

Eq. (A3) 

Cy 

V/£^ = 1 (box), 7r/6 (sphere) 

Eq. (A4) 

Cm 

^lo/^nt = 1 + 3/Ad^; 1 for highly supersonic flow 

Eq. (A5) 

C(^ 

Numerical coefficient in adopt c$ = l/(27r) 

Eq. (Al) 

M<b 

Magnetic critical mass 

Eq. (Al) 

Mb 

Alternate form of magnetic critical mass, M^j 

Eq. (AlO) 

M 

Sonic Mach number, Urms/cs 

Sec. 4.1 

Ma 

Alfven Mach number, Urms/'^^A 

Sec. 4.1 

nn 

Mean density of H nuclei 

Eq. (3) 

los 

Mean density of H nuclei in central Gaussian beam 

Eq. (21) 

-^H, beam 

Golumn density of H nuclei in central Gaussian beam 

Eq (21) 


Mass-weighted Alfven velocity, Brms ^/(^'Trp)^/^ 

Sec. 4.1 

"^^rms 

Mass-weighted 3D nonthermal velocity dispersion, v^3crnt 

Sec. 4.1 

(x); {x}j 

Average value of x in 100-clump sample; also averaged over time 

App. 5 

a 

Exponent for density scaling of magnetic field, B oc fig 

Eq. (3) 

0 

Plasma^, = 

Eq. (11) 


(Normalized) mass-to-flux ratio, MjM^ = 2 'kG^I‘^YiI B 

Eq. (A2) 

M‘i>,los 

Central-beam line-of-sight mass-to-flux ratio, 27rG^/^Sbeam/Mos 

Eq. (28) 


(l//i<i>}~ = harmonic mean mass/flux ratio = average flux/mass ratio 

Eq. (29) 


Mass-weighted total ID velocity dispersion, (Cg -|- 

Eq. (A5) 


Mass-weighted ID nonthermal velocity dispersion 

Eq. (A5) 

^beam 

Surface density in central Gaussian beam 

Sec. 6 


“ Subscripts: beam = central Gaussian beam, h = harmonic mean, H = hydrogen nuclei, los = line of sight, M = mass, med = median, 
t = time-averaged, V = volume-averaged, fi = microgauss 


Section 8. Theoretical results on the comparison of magnetic 
and gravitational forces are summarized in Appendix A. Fi¬ 
nally, Section 9 summarizes the results of our analysis of 
the line-of-sight data and our detailed comparison between 
simulations and observations of magnetic fields in clumps in 
molecular clouds. 


2 THE OBSERVATIONAL RESULTS OF 

CRUTCHER ET AL. (2010) 

The most comprehensive analysis of magnetic fields in the 
ISM is that of Crutcher et al. (2010), who analyzed Zeeman 
data from HI observations by Heiles & Troland (2004), OH 
observations by Crutcher (1999) and Troland & Crutcher 
(2008), and CN observations by Falgarone et al. (2008). The 
OH data comprises 54 sources, 20 of which are detected at 
more than 2 sigma; the CN data is from 14 sources, 8 of 
which are detected at more than 2 sigma. Crutcher et al. 
(2010) used a Bayesian analysis to self-consistently include 
upper limits on the field strengths with a statistical signifi¬ 
cance less than 2 sigma. 

The Zeeman signal is proportional both to the integral 
of the field along the line of sight and to the column den¬ 
sity of the molecule being observed. Let the line of sight be 


along the 2 direction. We assume that the molecule has a 
fixed abundance, so that the column density of the molecule 
is proportional to the mass column density. Then the mag¬ 
nitude of the Zeeman field measured by a given telescope 
beam is 


B\oS - T 7- 

/ pB,dV 

1 

/ BzdM 

pV 

/ 

^ beam 

^heam 

^ beam 


where Mbeam is the mass inside the beam and Bz is the com¬ 
ponent of the magnetic field along the line of sight. (Defi¬ 
nitions are collected together in Table 1.) Note that is 
a component of a vector and can be positive or negative, 
whereas values of the field that are not explicitly compo¬ 
nents of a vector, such as Bios, are field magnitudes. If the 
molecule being observed is excited only above some critical 
density, then Mbeam is the mass in the beam above that 
density. It is important to note that if the field is tangled, 
then Bios can be less than the average magnitude of B^. Of 
course, what is of greater interest is the value of the field in 
the absence of projection effects. 


Btot — 


A/beam 


BdM, 


( 2 ) 


and this is the focus of Crutcher et al.’s analysis. 

Crutcher et al. (2010) carried out a Bayesian analysis 
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to infer the distribution of Btot as a function of the mean 
density of hydrogen nuclei, hn, which extends over the range 
10 cm“®< hu < 10^ cm“® in their sample. They considered 
several different distributions and found that the best-fitting 
to their data is a uniform distribution of Btot from /Btot, max 
to Btot, max with 



r Bg (riH < no). 


-Stot, max — 

1 

(3) 


They found Bo = 10 /rG, no = 300 cm a = 0.65, and 
/ = 0.03 ~ 0. 

The density no dividing the constant density regime 
from the power-law one can be understood as the critical 
density, nu. crit, at which the thermal pressure in the clump 
balances the turbulent pressure of the ambient gas, 

Tturb = P^nt ~ ^totnu, crit^nT, (4) 

where p is the mean density, Xtot = ntot/hn is the total 
number density relative to hydrogen, and fee is the Boltz¬ 
mann constant. Since the magnetic field is about the same 
inside and outside the clump, it has no effect on the pressure 
balance. The observed clumps are in a variety of different en¬ 
vironments, and the minimum ambient pressure is that of 
the ISM. From Boulares & Cox (1990), the turbulent pres¬ 
sure in interstellar gas is Pturh/k « 1 X lO"* K cm Using 
T « 50 K and Xtot = 1.1 — 1 for atomic gas including he¬ 
lium, we find nn.crit ~ 200 cm“^, very close to the observed 
power-law relation starting point of 300 cm“^. 

All the data below no is from HI observations; most 
of the data above that is from OH and CN observations. 
Thus, the molecular data is inferred to have a power-law 
dependence on the density, 

Btot. max = 98hH®4 pG (molecular gas), (5) 

where hH.4 = hH/(10'* cm“®). It should be noted that the 
power-law dependence of Btot on hn is not apparent in the 
OH data alone (10® cm“® < hn 2 x 10'* cm“®) or the CN 
data alone (2 x 10® cm“® < hu 4, x 10® cm“®). A least 
squares fit of the OH data gives a slope of Bios versus hn 
of 0.40 ± 0.45 and that for CN data is 0.04 ± 0.67; that is, 
both the low density and high density molecular lines are 
consistent with either no density dependence or with the 
density dependence inferred from the combined data set. 
As they point out, a ~ the value expected for the con¬ 
traction of a spherical cloud with a weak, uniform magnetic 
field (Mestel 1966). An important implication of this model 
is that field strengths extend to very low values; for exam¬ 
ple, about 10 per cent of the clumps are inferred to have 
helds that are less than about 10 per cent of the maximum 
value expected at that density. The existence of very weak 
fields has been advocated as enabling the formation of discs 
around protostars (Krumholz et al. 2013). It should be noted 
that Btot, max is not an absolute upper bound on the field 
strength; both W3(OH) observed in CN (Falgarone et al. 
2008) and Mon 16 (Troland & Crutcher 2008) have fields 
exceeding Btot, max, although the latter result is uncertain 
since the field measurement is less than 2 sigma. 


3 NEW RESULTS FROM OLD DATA: 

ANALYSIS OF THE LINE-OF-SIGHT FIELD 

The results of Crutcher et al. (2010) show that the distri¬ 
bution of Btot/hg, with a ~ 0.65, is independent of density 
for molecular gas. In other words, while they presented their 
results as the scaling of the upper envelope of the total field, 
their best-fit distribution implies that the scaling applies to 
all the helds. We can test this directly: Since the line-of-sight 
held is a projection of the total held onto the line of sight, 
it follows mathematically that the distribution of Bios/hj), 
with a ~ 0.65, should also be independent of density in 
molecular gas. Indeed, a ht to the data on the line-of-sight 
helds in molecular gas from Crutcher et al. (2010), shown 
by the open green circles in Fig. 1, gives 

a(Bio.) = 0.64 ±0.13, (6) 

consistent with this expectation. Since this is indistinguish¬ 
able from the result of Crutcher et al. (2010), we shall adopt 
a = 0.65 as the observed value. We conclude that on aver¬ 
age, helds in molecular gas scale approximately as hn®®. We 
emphasize that the uncertainty of all the curve htting re¬ 
sults in this paper is at the 95 per cent conhdence level, i.e., 
about 2a. We choose to use a 95 per cent conhdence level 
since then any disagreement is likely to be a true difference 
instead of a statistical huctuation. We estimate the 95 per 
cent conhdence limits on the result of Crutcher et al. (2010) 
for a from their Figure 4: The probability that a < 0.59 or 
a > 0.75 are each less than 2.5 per cent, so their result is 
a — 0.65_Q gg. 

In Appendix E we discuss the effects of noise on the 
inferred value of a (Bios). There we show that since the noise 
in the low-density data (OH) is greater than that in the 
high-density data (CN), the value of a(Bios) inferred from 
the data is less than the true value by about 0.02. If the 
held is distributed isotopically with respect to the line of 
sight, then q;(Bios) = a. We conclude that the true value 
of a corresponding to the observed value a(Bios) = 0.64 is 
a — 0.66. This shift is small compared to the uncertainty, so 
it is not signihcant, but if the noise in future larger samples 
is asymmetric as in the present case this effect should be 
taken into account. 

Since the held scales as a power of the density, we dehne 


where B^ = B/(l pG). The quantity b for the line-of-sight 
held is denoted &ios and that for the total held btot- Hence¬ 
forth, we shall set a = 0.65 in this expression. For a sample 
of sources, it follows that 

(B^) = (b<®|) ~ {b)(h&®4®), (8) 

where the last step follows if B indeed scales as hn®®, so that 
b is not correlated with hn. We hnd that the average value 
of normalized line-of-sight held is {Bios, = 22.4, and 

that (Bios, 4 t)/(hH® 4 ) = 20.2. We conclude that the average 
line-of-sight held in the observed sample is related to the 
density by 

“ 'Values of o > 0.76 were not tested (Crutcher 2015, private 
communication), so our estimate of the upper confidence level is 
based on an extrapolation of the results for a Si 0.76 
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Figure 1. Dependence of the observed field strength, -Bios, 
density, nn. The convolved, mass-weighted, line-of-sight magnetic 
field strength, Bios, along the three cardinal axes vs volume den¬ 
sity in the central-beam regions of clumps from the strong-field 
simulation (black pluses) are plotted with the observed HI (green 
diamonds) and OH-hCN sources (green circles) from Crutcher et 
al. (2010). The dashed line is the best power-law fit to the sim¬ 
ulation data. The large dispersion is due primarily to projection 
effects. 

fiG. (9) 

As noted after Equation (B5) in Appendix B, the average 
total field is just twice the average line-of-sight field. We 
therefore conclude that the average magnetic field in the 
molecular clumps in the study of Crutcher et al. (2010) is 

{Btot(hH)) « 424®! pG. (10) 

Insofar as the clumps in this study are broadly representa¬ 
tive, this is the average field molecular clumps in the inter¬ 
stellar medium. 

The Bayesian analysis of Crutcher et al. (2010) en¬ 
ables one to infer the range of possible 3D magnetic fields 
that are consistent with the observed distribution of line-of- 
sight fields. However, it is also possible to infer the prop¬ 
erties of the fields directly from the line-of-sight data. As 
shown in Appendix B, the mean value of the total field 
can be inferred from the mean value of the line-of-sight 
field, (Btot) = 2(Hios), and the density dependence of the 
field can be directly measured as well. Furthermore, instead 
of comparing the inferred 3D fields with a model distribu¬ 
tion for those fields, we shall convert the 3D distribution 
to the line-of-sight distribution and compare that directly 
with the data, thereby eliminating any uncertainties intro¬ 
duced by the statistical inference of the 3D properties of 
the field. Since the physical conditions we have chosen for 
our simulation are appropriate for molecular clouds, we shall 
compare our results only with the 68 molecular clumps in 
the Crutcher et al. (2010) sample. Our sample therefore 
excludes three HI sources with densities in the power-law 
region (hu > 300 cm“®) found by Crutcher et al. (2010); 
in fact, these sources have densities ^ 5000 cm“®, overlap¬ 
ping those of molecular clouds. Only 28 of the 68 molecular 
clumps in this sample have magnetic fields detected at a sig¬ 



Figure 2. Model fits to the distribution of Bios/n-H®® for the 
observed OH-I-CN sources: lognormal (blue solid line), exponen¬ 
tial (black dot-dash line), and flat distribution (red dashed line). 
The differences among the three models are small (see Section 3 
for discussion). More observations are required to determine the 
actual distribution. 

nificance greater than 2 sigma. For the measurements that 
are less than 2 sigma, we use the measured value, which 
is appropriate provided that the errors are statistical and 
not systematic. Since this is also a necessary condition for 
the validity of the analysis of Crutcher et al. (2010), noth¬ 
ing is lost by this procedure. This method may be thought 
of as analogous to stacking optical or X-ray data to bring 
out results that are buried in the noise in individual frames. 
However, there is a difference in the present case: whereas 
photon noise can result in the observed value being less than 
the noise, it is always positive; for Zeeman observations, on 
the other hand, measurement errors could change the sign 
of the observed magnetic field if the noise is larger than the 
true field strength. In Appendix C we discuss a numerical 
experiment that demonstrates that this effect is small. 

We have considered three theoretical models for the 
distribution of field strengths: (1) a uniform distribution 
of total normalized field strengths from /Hmax/uH^I to 
Bmax/hH ®4 (Crutcher et al. 2010); (2) a log-normal distribu¬ 
tion of total normalized field strengths; and (3) an exponen¬ 
tial distribution of line-of-sight normalized field strengths. 
Explicit expressions for the total fields and line-of-sight fields 
for each are given in Appendix B. The distribution of the ob¬ 
served values of Hioa, /j/h-H®! are compared with these models 
in Fig. 2. The striking feature of this figure is that the three 
very different distributions of total field strength, Btot, look 
remarkably alike when viewed along the line of sight, provid¬ 
ing graphic evidence of the difficulty in accurately inferring 
the distribution of Btot from observations of the line-of-sight 
field. 

We conclude that the log normal and exponential mod¬ 
els provide slightly better fits to the data than the uni¬ 
form distribution. This does not contradict the conclusion 
of Crutcher et al. (2010): For the molecular data, they con¬ 
sidered three cases: (1) OH data alone with B independent 
of density; they found that a uniform distribution was bet¬ 
ter than one obtained from a numerical simulation, which 
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Figure 3. The best-fitting power-law relation between 
and the telescope-beam size rteam of 34 OH and and 14 CN 
sources (open circles) from Troland & Crutcher (2008) and Fal- 
garone et al. (2008). The three solid circular symbols are from the 
100-clump strong-field sample, using 3 different beam sizes (0.025 
pc, 0.0125 pc, and 0.00625 pc) at t = OAtg. See the discussion in 
Sections 2 and 6.1. 


was approximately log normal; (2) OH -|- CN with an as¬ 
sumed density dependence B oc here they found that 
a uniform distribution was better than a delta function, but 
they did not consider a log normal; and finally, (3) all the 
data, including HI, with a density dependence that was de¬ 
termined from the data, B oc they did not consider a 

log normal fit to the data in their paper. Thus in no case 
did they compare a log normal vs. a uniform distribution 
for the molecular data with B oc as we have done here. 
More accurate Zeeman measurements and more sources are 
required for a better comparison. 

Before concluding our discussion of the observations, 
we note that we have found evidence that the measured 
values of Sios/h-H®* for the OH sources studied by Troland 
& Crutcher (2008) and the CN sources studied by Falgar- 
one et al. (2008) depend on size of the beam, rbeam, used 
to study them. (We do not know the beam sizes for the 
sources listed by Crutcher (1999), which came from a vari¬ 
ety of references.) In Fig. 3, we plot the measured values of 
Bios/fin®* for these 48 sources as a function of rbeam. The 
correlation coefficient is 0.271 and the probability that there 
is no correlation isp = 0.062. The correlation is weak, as the 
scattered data points show in Fig. 3, but p is small enough 
that this correlation is marginally significant. Further obser¬ 
vations are needed to determine if this correlation is real. 


4 SIMULATIONS 

We have carried out two large-scale AMR simulations to 
study the formation and structure of Infrared Dark Clouds 
(IRDCs) (Li et al. 2015, in preparation). These simulations 
use the ORION2 adaptive mesh refinement (AMR) code (Li 
et al. 2012a) with ideal MHD and self-gravity to follow the 
evolution of driven turbulence on a base grid of 512® with 
two levels of refinement, with each level having a factor of 


2 increase in resolution. We adopt periodic boundary con¬ 
ditions. We drive the system for two crossing times with 
no self-gravity in order to produce self-consistent turbulent 
density and velocity fields. The first crossing time is at the 
base level and the second has two levels of refinement. The 
refinement thresholds for pressure jumps (thermal pressure 
-I- magnetic pressure) and shear flows are chosen so that 
about 12 — 15 per cent of the volume is refined at all times 
during the second crossing time. This ensures that before 
gravity is turned on, the high-mass end of the probability 
density function (PDF) of volume density is equivalent to 
that which would be obtained in a 2048® single-grid simula¬ 
tion (Li et al. 2012a). As we shall see in Section 4.2 below, 
the simulation box has a size of 4.55 pc for fiducial values 
of the parameters, so the cell size at the finest resolution is 
2.2 X 10“® pc. The driving is at the largest scale, with wave 
numbers fc = 1 — 2, and it uses the Mac Low (1999) recipe. 
After two crossing times, we turn on gravity and set t = 0 
at this moment. The driving continues with a constant en¬ 
ergy injection rate that maintains the system globally at a 
constant sonic Mach number, M. Dense regions inside the 
system undergo gravitational collapse, so it is necessary to 
add the Jeans condition for refinement (Truelove et al. 1997); 
we adopt a Jeans number of 1/8. The finest resolution re¬ 
mains unchanged from the first part of the simulation, and 
no sink particles are introduced. 


4.1 Scaling Relations for a Turbnlent Box 


Simulations of isothermal MHD turbulence are scale-free 
and are defined by two dimensionless parameters, the three- 
dimensional (3D) thermal Mach number, M = Urms/cs, 
where Vrms is the mass-weighted rms velocity, and the 
Alfven Mach number, Ma = v^ma/vA, where va = 
Brms.v/(4:7rp)®fo is the mass-weighted Alfven velocity and 
Brms.v = is the square root of the volume- 

weighted mean-square magnetic field. Alternatively, the 
strength of the magnetic field can be characterized by the 
plasma-/? parameter. 


SnpCa 

B® 

■^rms, V 


= 2 



( 11 ) 


In order to set the scale, we need three dimensional 
relations (McKee et al. 2010). Interstellar molecular gas 
is observed to have a temperature T = lOTi K with 
Ti ~ 1, corresponding to an isothermal sound speed Cs = 
0.188r/^^ km s“®; this gives one relation. Most interstellar 
gas obeys the turbulent linewidth-size relation (McKee & 
Ostriker 2007), 


ant = 0.72aCRli^ = O.blalX^l^ km s"®, (12) 


where ant is the mass-weighted ID nonthermal velocity dis¬ 
persion in a sphere of radius R, which we take to be the 
same as in a box of size io = 2B, and where a*;, ~ 1 al¬ 
lows for deviations from the typical linewidth-size relation. 
When self-gravity is included, a third dimensional param¬ 
eter enters, the gravitational constant, G, together with a 
corresponding dimensionless parameter that measures the 
effect of self-gravity, the virial parameter, Ovir — 2x kinetic 
energy/(gravitational energy) (see Equation (A6)). We as¬ 
sume that the turbulent box is large enough that the gas is 
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Table 2. System Properties for the Strong- and Weak-Field Simulations 




Strong 



Weak 


Initial" 

t = 0® 

Final" 

Initial" 

t = 0® 

Final" 

Alfven Mach number, A^A(^rms) 

1 

0.85 

0.83 

10 

1.86 

1.61 

Plasma /3 = Sirpoc^ / 

0.02 

0.015 

0.014 

2.0 

0.069 

0.055 

Mass-to-flux ratio, z 

1.62 

1.62 

1.62 

16.2 

16.2 

16.2 

Avg magnetic field (/^G), {|B^|}y 

31.6 

36.2 

36.8 

3.16 

14.9 

16.2 

RMS magnetic field (^G), Brms, 

31.6 

37.2 

38.3 

3.16 

17.1 

19.2 

D /-2/3 

^rms, 4 

151 

177.4 

182.7 

15.1 

81.6 

93.0 


^ The initial field is uniform so that Bq = (|i?|)y = firms then. 

Gravity is turned on at t = 0, after the system is driven for two crossing times. 

The final values are time averages over (0.4 — 0.64)tfi for the strong-field case and (0.5 — 0.94)f£f for the weak-field case. 


highly supersonic so that the mass-weighted ID total veloc¬ 
ity dispersion, ctid, obeys nfjj = cr^t +Cs ~ For a highly 
supersonic, fully molecular gas, the physical parameters of 
the turbulent system—the size of the turbulent box, £o, the 
flow time (or crossing time), tf, the mass of the box. Mo, 
and the column density, A^h = nnto —can then be expressed 
as (see the Appendix in McKee et al. 2010): 


to 

if 


2 M'^c? f MMi 

= = 0.0455 I I pc. 


3 (cr^c/ Ipc) 


= 2.36 X 10® ' 


1/2 


yr, 


riH 

= 9.6 X 10“ 

Mo 

= 0.311 f- 


v«- 

A^h 

= 1.34 X lO' 


4 \ 


rtvirM^Tl ) 

> M«, 


22 I '^pc 
rtvir 


(13) 

(14) 

(15) 

(16) 
(17) 


It should be noted that the only assumptions entering the 
derivation of these relations is that the Mach number is large 
and that the velocity dispersion scales as ■ Observe that 
the size scale for a simulation of self-gravitating (civir ~ 1) 
interstellar molecular gas (Ti ~ 1) that obeys the standard 
linewidth-size relation (o-pc ~ 1) is set solely by the Mach 
number, A4. Equation (17) shows that the column density of 
such gas is independent of the Mach number and is therefore 
approximately constant, one of the three relations discovered 
by Larson (1981). With different coefficients that allow for 
the difference in shape between an interstellar cloud and a 
turbulent box, these relations apply to actual interstellar 
clouds as well. 

The magnetic field is given by 


(18) 

where hn.s = hH/(1000 cm~®). Magnetic fields in gas that 
obeys the standard linewidth-size relation can be signfi- 
cantly less than 30 fiG only if self-gravity is weak (ovir »1) 
and/or if the motions in the gas are highly super-Alfvenic 
{Ma > 1). 

The relative importance of magnetic fields and grav¬ 
ity is discussed in Appendix A. Clouds with masses exceed¬ 


ing the magnetic critical mass, M^, can undergo gravita¬ 
tional collapse; those with smaller masses cannot. We adopt 
M<s. =(magnetic fiux)/27rG^^^, which is exact for a sheet 
threaded by a perpendicular field (Nakano & Nakamura 
1978). The mass-to-fiux ratio normalized to the critical value 
is then 


— ^ r, ^1/2^1 


3.80 / A^h \ 

Bfj, V 10^1 cm“^ / 


(19) 


4.2 Model Parameters 

The models we consider here have A4 = 10, T = 10 K, and 
rtvir = rtp,, = 1, so that the size of the simulated turbulent 
region is to = 4.55 pc, the number density hn = 960 cm“®, 
the total mass Mq = 3110 Mq, and the initial column den¬ 
sity of the system A^h = 1.34 x 10^^ cm“^. The free-fall time 
of the system based on the mean density, 

/ o \ 1/2 

= ( ^ ) = yr, (20) 

is 1.4 X 10® yr, which is 0.59tf. 

We present two simulations in this paper, one with a 
strong initial mean field (A4 a.o = 1, Ms.o = 1-62), which 
is only slightly supercritical, and one with a weak initial 
mean field (A4 a,o = 10, = 16.2), which is initially 

very supercritical. The parameters of the runs are given in 
Table 2. The initial magnetic field is along the 2 -axis. We ran 
the strong-field simulation to 0.64tff « 900,000 yr , which 
provided enough data to time-average the parameters that 
do not vary systematically with time for comparison with 
observation. We ran the weak-field model somewhat longer, 
to t ~ 0.94tff, since the maximum column density in our 
simulation is comparable to that in the observed sample at 
that time. A detailed discussion of the time evolution of the 
whole system will be given in a subsequent paper (Li et al. 
2015, in preparation); here we focus on the properties of the 
magnetic fields in dense clumps during the self-gravitating 
phase of the simulations. 
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5 THE 100-CLUMP SAMPLE FOR 

COMPARISON WITH OBSERVATION 

Since the magnetic fields in the Crutcher molecular clump 
sample are measured in condensations in molecular clouds, 
we select a sample of dense clumps from our simulation. The 
most massive dense clumps are the most readily observable, 
so we select them; in order to have a good statistical sample, 
we choose the 100 most massive dense clumps and term this 
the TOO-clump sample.’ It Is important to note that during 
the course of the simulation, clumps accrete matter and can 
merge, so the clumps in this sample at one time are not 
the same as those at a different time. We identify the dense 
clumps from the simulations using a modified CLUMPFIND 
algorithm (Williams, De Gens, & Blitz 1994). The clumps 
are defined as connected regions with a density greater than 
the mass-weighted median density of the turbulent region, 
(h-H)M, med; half the mass in the region has a higher density 
than this and half a lower one. No cell with density below 
this threshold will be included as part of the clumps that 
are identified. This density rises from about Snn at t = 0 
to about 6nH at the end of the simulations as a result of 
gravitational contraction, although there is no requirement 
that the clumps be gravitationally bound. We note that the 
clumps for which Zeeman observations are available are not 
necessarily bound either (e.g., Falgarone et al. 2008). 

To assemble a sample of simulated clumps for compari¬ 
son with observations, we impose requirements on the reso¬ 
lution, density, and mass as described in Appendix C. Reso¬ 
lution: The minimum radius of the clumps must exceed the 
minimum radius of the OH and CN sources, which is 0.03 pc; 
all clumps closer than 0.06 pc are merged together. Density: 
As noted in Section 2 above, the power-law behavior of the 
field begins at a density no — 300 cm“® that can be under¬ 
stood as the critical density at which the thermal pressure 
of the gas balances the typical turbulent pressure of the in¬ 
terstellar medium. Our simulation is meant to represent the 
interior of a GMC, so the pressure is significantly higher; 
for a molecular gas at T = 10 K, the corresponding criti¬ 
cal density is ~ 3.2 x 10"^ cm“®. We require the clumps to 
have densities slightly higher than this, nn > 4 x 10"^ cm“^. 
With this choice of density threshold, the simulated clumps 
have densities closer to those of the CN sources than of the 
OH sources. Mass: As noted above, we select the 100 most 
massive clumps that satisfy the above two requirements. 

Although we have selected our clumps from simulations 
based on the typical conditions observed in molecular clouds 
(in particular, the simulations have A4a,o = ctvir = 1), there 
are important differences between our simulations and re¬ 
ality. The molecular clumps in the Crutcher et al. (2010) 
sample are from a variety of clouds with varying line width- 
size relations, virial parameters, and mass-to-flux ratios. The 
initial conditions of our two simulation models have only one 
line width-size relation and one virial parameter. We know 
exactly the duration of the gravitational collapse phase of 
the simulations, but the cloud clumps observed can be of 
any age. Our simulations are isothermal, but there are tem¬ 
perature variations in real clouds. Observed cloud clumps 
are obtained from OH and CN emission-line observations, 
but our clumps are identified using CLUMPFIND and the 
rules listed above. Some observed clouds have masses of the 
same order as the total mass of our whole simulated region. 


whereas most of the clumps in the simulations are less than 
IOMq. On the other hand, the densities in our clump sam¬ 
ples are comparable to to those in some OH sources and in 
all the observed CN sources. Despite these differences, we 
shall show in the following sections that there is excellent 
agreement between the properties of clumps in the simula¬ 
tions and those of the observed clumps. 

To compare with observations or to compare the two 
models, we use the mean, denoted (x) for a quantity x, 
and median, (a;)med, from the 100-clump samples. All ob¬ 
served magnetic fields are mass weighted. When ambiguity 
is possible, mass-weighted quantities are denoted by {x}j^ 
and volume-weighted ones by {x}y. When comparing time- 
dependent quantities, we use data from t = 0.57tff for the 
strong mean-field simulation and 0.94tfr for the weak mean- 
field simulation since the maximum column densities in the 
100-clump samples are comparable to the maximum value in 
the Crutcher et al. (2010) sample at those times. When com¬ 
paring time-independent quantities, we evaluate the time 
average, denoted by the subscript t. For example, the time- 
averaged mean is {x}t- For the strong mean-field model, we 
determine the time averaged value from data at t = 0.4, 0.5, 
0.57, and 0.64tff; for the weak mean-field model, which takes 
somewhat longer to reach an approximate steady state, we 
use data at t — 0.5, 0.6, 0.7, 0.8, 0.9, and 0.94t£f for the time 
average. 

We did not use sink particles in these simulations, so 
some of the clumps may violate the Jeans condition, which 
is that the Jeans length must be resolved by enough cells to 
prevent artificial fragmentation (Truelove et al. 1997). In our 
100-clump sample, we require the clumps to be at least 0.03 
pc in radius, which eliminates the possibility of separating 
small artificial fragments into separate clumps. To quantify 
how much clumps that violate the Jeans condition could af¬ 
fect our results, we removed the clumps in which cells with 
more than 25 per cent of the mass violated the Jeans condi¬ 
tion and replaced them with clumps of lower mass so as to 
maintain a 100-clump sample. We found that the power-law 
index describing the relation between B and nn discussed in 
the next section changed by only 0.02, well within the un¬ 
certainties. It thus appears that our results are not sensitive 
to artificial fragmentation. 

In Appendix D we show that the results of our simu¬ 
lation for the 100-clump strong mean-field sample are con¬ 
verged by comparing with results of a simulation with 3 
levels of refinement. The convergence test was carried out 
at t = O.dtff, the latest time that we could afford to run 
with 3 levels of refinement. Almost all of the quantities from 
the 2-level and 3-level lOO-Clump samples in Tables A.l and 
A. 2 are converged within the error of the means. The one 
exception is the mass-weighted rms field, Srms.M, which has 
a difference slightly larger than the error of the individual 
means, although within the rms error. We repeated the test 
described above, in which we analyzed the results after elim¬ 
inating the clumps that violated the Jeans condition at the 
25 per cent level, for both the 2-level and 3-level 100-clump 
samples and found that the results are still converged, show¬ 
ing once again that even though there are cells that violate 
the Jeans condition, this has no significant effect on our re¬ 
sults. 
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Figure 4. (a) The CN image of molecular cloud W3(OH) using CARMA (Hakobian & Crutcher 2011), showing clumpy structure, (b) 
The column density map of a massive clump from the strong-field simulation at 0.57%, showing structure inside the clump. 


6 COMPARISON OF THE STRONG 

MEAN-FIELD MODEL WITH 

OBSERVATIONS 

We wish to compare our simulations both with each other 
and with observation. In Section 8 we shall see that although 
many of the results from the two simulations are similar, the 
weak-held simulation does not reproduce the observed cor¬ 
relation of the magnetic field orientation on large and small 
scales (Li et al. 2009). Since the strong-if eld model is favored 
by observation, in this section we shall compare the results 
of that simulation with the observations of Crutcher et al. 
(2010). Equations (1) and (2) show that we need the mass- 
weighted fields averaged over the telescope beam in order to 
compare with the observations. Zeeman observations require 
long observing times, so that only the single region of a given 
clump with the strongest emission in the Zeeman molecule is 
observed (Crutcher, private communication). Fig. 4a shows 
a CN map of a clump in the molecular cloud W3(OH) ob¬ 
served by Hakobian & Crutcher (2011). Fig. 4b is a column- 
density map of one of our massive clumps. Like the clump 
in W3(OH), it shows considerable internal structure. The 
complexity of the structure in both the observed clump and 
the simulated one in Fig. 4 shows that, like the observers, 
we are obtaining only the average of one component of the 
field in the clumps. 

As discussed in the selection of the 100-clump sample in 
Section 5, the smallest known beam size of the CN and OH 
sources in Crutcher et al. (2010) is about 0.025 pc in radius. 
We ‘observe’ our clumps also using a Gaussian beam with 
radius of 0.025 pc. For column density, we create column 
density maps of a clump along each cardinal axis {x, y, and 
2 -axes) and convolve the maps using a Gaussian beam of 
radius 0.025 pc. The peak column density of each convolved 
map is the ‘observed’ column density of the clump along 
that direction: Ah, beam is the number of H nuclei per unit 
area and Ebeam = (2.34 x lO”^'^ g)A h, beam is the mass per 


unit area in the beam. To obtain the observed magnetic 
field. Bios, we create a density-weighted map of the line-of- 
sight magnetic field integrated along each cardinal axis first, 
convolve the map with the same Gaussian beam, and then 
take the magnitude. The ‘observed’ field strength is the value 
of Rios at the location of the peak ‘observed’ column density 
in the convolved column density map. We thus obtain 300 
measurements of Ah, beam and the corresponding Bios from 
our 100-clump sample. The mean volume density of each 
clump is obtained from the measured column density as 

Uh, los — Ah, beam/2rclump, (^1) 

where rdump is the mean radius obtained from the contour 
of the half maximum of Ah. 


6.1 Density Scaling of the Magnetic Field 

One of the primary results obtained by Crutcher et al. (2010) 
is that in molecular gas there is a power-law relation between 
the value of Rtot, max that they infer (but do not directly 
measure) and the density, Rtot.max oc hg, with a = 0.65. 
As we discussed in Section 3, this relation applies to the ob¬ 
served line-of-sight fields and to the total fields as well. In 
many cases the density of the molecular gas was estimated 
as in Equation (21), although in some cases the density was 
inferred from observations of density-sensitive lines. At low 
densities, hn.ios 300 cm“®, they found that the clumps 
(which are all HI at these densities) form a horizontal branch 
with low field strengths, Rios ^ 10 yG, that do not scale with 
density. In creating our 100-clump sample, we eliminated the 
clumps in the horizontal branch since our focus is on clumps 
in molecular clouds, which are observed to be in the power- 
law regime. The density range of the central-beam region of 
the simulated clumps is in the range (5 x 10"* — 1.1 x 10^) 
cm“®, which is somewhat above the density range of the OH 
observations but covers the range of the GN observations. 
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Figure 5. Power-law relation between total field strength, Btot, 
and volume density of the central beam region in the strong-field 
100-clump sample at 0.57ffi (open circles and solid line). The 
power-law index relating Btot to njj at that time is oo.STtff = 
0.62 it 0.11. The time-averaged [(0.4 — 0.64)%] value from the 
strong-field 100-clump samples is shown by the dashed line la¬ 
beled {a)t = 0.70 it 0.06. Clumps that have column densities of 
A^h ^ 2 X 10^^ cm“^ at t = 0, prior to the time that self-gravity 
has acted, are shown by pluses. At that time, the magnetic field 
strength of clumps in the simulation is independent of the vol¬ 
ume density, similar to the horizontal branch (HI cloud clumps) 
in the Btot vs njj diagram in Fig. 1. The dash-dot line indicates 
the average field strength of these clumps. The power-law branch 
of Btot vs. nj] is the result of gravitational collapse of clumps to 
higher density and magnetic field strength. 


As discussed in Section 5 above, the transition from the hor¬ 
izontal Bios(hH) relation to the power-law regime begins at 
a much higher density in our simulations than in the ob¬ 
servations because the transition occurs when the thermal 
pressure in the clumps is comparable to the ambient turbu¬ 
lent pressure, and this occurs at a much higher density in 
the molecular clouds we are simulating than in the diffuse 
ISM that is included in the observations. 

We now demonstrate that our simulations show that 
Btot oc hg, with a value of a consistent with that found 
by Crutcher et al. (2010) and obtained in Section 3. We 
begin with a discussion of how we obtain Btot from our sim¬ 
ulations. Since each core in the simulation gives us three 
different line-of-sight fields to analyze, we infer Btot based 
on observing the field in the three cardinal directions, 

Btot — (Bios, 3:, B\oB,y, Bios, 2 ). 

Thus, Btot in Equation (2) is the field associated with the 
mass in the actual telescope beam, which we have assumed 
is on the z axis, whereas Btot in Equation (22) is the field 
that would be measured by an Earthbound observer and two 
hypothetical observers along the x and y axes, who would 
observe slightly different volumes of gas unless the emission 
were unresolved. We have shown that the distributions of 
Btot = I Btot I calculated with Equations (2) and (22) are 
very nearly the same. Furthermore, the time-averaged me¬ 
dian value of Btot/Bios using this procedure is 2.05, very 
close to the expected value of 2 for the case in which the 


same region of the clump is used to determine the relation 
between Btot and Bios. 

The values of Bios at t = 0.57% in the strong mean-field 
simulation are plotted as a function of hn in Fig. 1 along 
with the line-of-sight field measurements of Crutcher et al. 
(2010). The distribution of the line-of-sight fields versus the 
mean density determined along each line of sight in the 100- 
clump sample fully covers the distribution of observed CN 
clumps in terms of hn and the top envelope of the distribu¬ 
tion is consistent with a power-law relation between Bios and 
hn, similarly to the power law obtained from the Bayesian 
analysis of Crutcher et al. (2010). The distribution is also 
widely scattered, similar to that of the observed CN and 
OH clumps. Below, we will see that (Bios/hH^*)med for the 
100-clump sample is close to the observed value. 

The relation between the field strength and the den¬ 
sity is much clearer for Btot since the scatter due to pro¬ 
jection along the line of sight is removed. For the density, 
we take the mean value in the central beam region along 
the three cardinal axes. As shown in Fig. 5, the strong-field 
100-clump sample at t = 0.57tfi shows a clear power-law 
relation between Btot and hn with a = 0.62 ± 0.11, con¬ 
sistent with the value inferred from observation. This result 
is explored further in Fig. 6, which shows a at several dif¬ 
ferent times (results of the weak mean-field model are also 
shown). Note that a is approximately independent of time 
once self-gravity has had time to act; by contrast, a can¬ 
not be determined at t = 0 since Btot and hn are uncor¬ 
related then, as shown in Fig. 5. The time-averaged value 
of a over the period (0.4 — 0.64)tff is (a)t = 0.70 ± 0.06, 
consistent with the value of a = 0.65 from Crutcher et al. 
(2010). The line-of-sight fields in the simulations exhibit a 
power-law relation with density as they should, since they 
are the projection of fields that obey such a relation; for the 
strong-held case, we find (a(Bios))t = 0.61 ±0.07, consistent 
with the value a (Bios) = 0.64 from a least-squares fit to the 
data of Crutcher et al. (2010). The results for the power-law 
dependence of the magnetic field on the density are sum¬ 
marized in Table 3. The important conclusion to be drawn 
from these results is that the observations and simulations 
are consistent, and both show that Btot and Bios scale as 
within the errors. 

In addition to determining the power law in the relation 
between the field and the density, it is important to deter¬ 
mine the coefficient. In Appendix A we show that Bfri^^ 
is proportional to a magnetic critical mass, Mb, that does 
not depend explicitly on the size of the clump (Eq. AlO). 
Equivalently, for gravitationally bound (cuvir ~ 1), super¬ 
sonically turbulent (A4 ^ 1) clumps, B/n^^ oc 
is inversely proportional to the mass-to-fiux ratio y,s>. These 
results change very little when the exponent 2/3 is replaced 
by the observed value of a = 0.65. Thus, the coefficient in 
the power-law relation depends on the mass-to-fiux ratio. 

We have collected the time-averaged mean and median 
values of both Btot, and Bios, y/nu^i from our sim¬ 
ulations together with the observed values in Table 3. Here 
we briefly compare the strong-field central-beam values with 
the observations; we defer discussion of the the whole-clump 
values to Section 7 and the weak-field values to Section 8. 
The time-averaged median value of Bios, is about 75 

per cent of the observed value. This agreement is extraor¬ 
dinary, both because of the many differences between the 
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Figure 6. Approximate time independence of key quantities: The power-law index a, normalized magnetic field strengths, fJtot/njr 
and Sios/fiH’ line-of-sight mass-to-flux ratio, from the 100-clump samples in the two simulations at different times in the 

gravitational collapse phase. The strong-held results are in the top row (a-c) at times 0.4, 0.5, 0.57, and O.Odtjf; the weak-held results are 
in the bottom row (d-f) at times 0.5, 0.6, 0.7, 0.8, 0.9 and 0.94 tff. Values at t = 0 are shown for Btot/n'^ and (a) The horizontal 

solid line is the time-averaged mean, {a)t, over this period of time, and the two dot-dashed lines show the uncertainty of (o)t. (b) The 
means (open circles) and medians (solid circles) of Stot/uH®® and those (squares) of Bios/h^®®. Comparison with the time-averaged 
means (solid lines) and medians (dashed lines) of and Sios/h^®® (including t = 0), respectively, show that they do not vary 

much from the initial value at t = 0. The value of is about half that of as expected from a random distribution 

of Held orientations, (c) The harmonic means (open circles) and medians (solid circles) of show a slow increase with time. The 

time-averaged harmonic means (solid lines) and medians (dashed lines), excluding i = 0, are shown for comparison, (d), (e), and (f) are 
the same as (a), (b), and (c), respectively, but for the weak-held model. 


simulated clumps and the observed clumps discussed in Sec¬ 
tion 5 and because the simulated clumps and the observed 
clumps each have a large dispersion in Bios (see Fig. 1). On 
the other hand, the average observed value of Bios/h-H®| is 
about twice that in the simulation. 

What causes the power-law relation between the mag¬ 
netic field and the density? In our simulation, it is clearly 
self-gravity. Fig. 5 shows the 100-clump sample from the 
strong mean-field simulation at t = 0.57tff in the hu- 
Btot plane, and the figure also shows all the clumps with 
A^h = 2 X 10^^ cm“® before gravity is turned on {t — 0). 
The power-law tail is present only after self-gravity has had 
time to act. The same is likely to be true in the obser¬ 
vations: The HI clumps are not self-gravitating and have 
Btot independent of density, whereas many of the molecular 
clumps, which do show Btot oc hn, are self-gravitating (e.g., 
Crutcher 1999). In Equation (4) we showed that the power- 
law region begins at about the point at which the turbulent 
and thermal pressures are the same; Krumholz & McKee 
(2005) showed that this corresponds to the condition for the 


clump to be self-gravitating, provided it is embedded in a 
self-gravitating cloud. The effect of self-gravity on the 100 
clumps can be assessed by determining if they are collaps¬ 
ing and, if so, how the collapse velocity compares with the 
free-fall velocity. We compare the mean value of the density- 
weighted infall velocity averaged over a sphere of radius r 
to the free-fall velocity there, Uff(r) = [Gm(r)/r]^^^, 

r pVrdA 

vs{r) f pdA' 

We have evaluated this ratio for the 100 clumps at O.bTta 
at a radius of 0.05 pc. The median value is 0.29, showing 
that these clumps are contracting at a significant fraction 
of the free-fall velocity and confirming that self-gravity is 
important. 

Crutcher et al. (2010) pointed out that a ~ | is consis¬ 
tent with the suggestion by Mestel (1966) that the spherical 
contraction of a cloud with a weak magnetic field would 
lead to B oc . In Section 7.2, we shall see that the me¬ 
dian mass-to-flux ratios of the 100-clump samples from the 
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Table 3. Observations Compared with Time-Averaged Properties of 100-Clump Samples 


Strong field“ Weak field" Observed" 


Central Beam^ 


(a(Bi, 

os)}t 

0.63 

± 

0.07 

-H 

00 

o 

0.05 

(Bios, 

/fjO.SSl 

m/'‘h,4A 

11.1 

± 

0.4 

9.2 ±C 

1.4 

(Bios, 


8.9 



7.1 


{a{Btot))t 

0.70 

± 

0.06 

0.57 ± 

0.05 

(Btot, 

/js0.65\, 

M/^H iV 

21.1 

± 

0.8 

18.2 ± 

1.6 

(Btot, 


18.7 



15.7 



0.64 ±0.13 
22.4 ±3.3 
11.9 

44.8 ± 6.6 


Whole Clump 


(a(Bi, 

os))t 

0.67 

± 

0.10 

0.65 

± 

0.08 

(Bios, 

/n“'®5V 

m/'‘h,4A 

10.7 

± 

0.4 

11.9 

± 

0.4 

(Bios, 

M/4'®|)med,t 

8.8 



8.6 



(o(Btot)}t 

0.70 

± 

0.06 

0.69 

± 

0.07 

(Btot, 

/js0.65\, 

M/^H, 4 A 

21.9 

± 

0.4 

23.7 

± 

0.8 

(Btot, 


19.1 



19.7 




" See Table 2. The time averages cover (0.4 — 0.64)tfif in the strong-field case and (0.5 — 0.94)tff in the weak-field case. 

Both i?ios and nn are convolved with a Gaussian beam 0.025 pc in radius (see Section 6). 

" Based on data from Crutcher et ah (2010). The time averages of the quantities in the table do not apply to the observed results. 

These 95 per cent confidence limits are inferred from Fig. 4 of Crutcher et ah (2010) under the assumption that the count above 
a = 0.76 is small. 


two simulations are around 2, with the strong-field model 
being somewhat less than 2 and the weak-field one some¬ 
what greater. The median magnetic field is therefore not 
very weak compared to gravity. However, the contraction of 
the clumps is approximately spherical: The time-averaged 
median value of the aspect ratio ZjR of the 100-clump sam¬ 
ples from the strong and weak mean-field models at the end 
of the simulations are 1.19 and 1.08, respectively. (Recall 
that Z is along and R is perpendicular to the mean field 
direction of the clump.) In fact, all clumps in both samples 
have 0.5 < Z/R < 2 at all times t ^ 0 and the medians re¬ 
main close to unity throughout the simulations, so that their 
collapse is approximately spherical. Note that the range of 
ZjR implies that there are clumps with prolate, oblate, and 
triaxial shapes in the 100-clump samples, as seen in other 
ideal MHD simulations with turbulence and self-gravity (e.g. 
Gammie et al. 2003; Li et al. 2004). 

Independent of the strength of the field, a necessary 
condition to obtain B oc is that the coefficient in the 
B{nH) relation must be about the same for all the clouds. 
The results of our simulation show that the ratio of the 
time-averaged dispersion of btot = Btot, which we 

label (j(&tot), to the median value is 


(btot, med)^ 

which is small enough so that this is approximately true. 
This ratio for bios is 0.60; it is larger because of the additional 
dispersion added by projection effects. 

We conclude that our simulations show that the density 
dependence of the field found by Crutcher et al. (2010) is 
due to gravitational contraction, that this contraction is ap¬ 
proximately spherical, and that different clumps have similar 
values of the ratio . These characteristics are all con¬ 

sistent with the formation of molecular clumps due to the 
gravitational contraction of clouds that initially had similar 


values of B and hu (e.g., in the cold neutral medium of the 
ISM) under the assumption that the field is so weak that 
it does not affect the contraction. However, although the 
field is weak enough that almost all the clumps are magnet¬ 
ically supercritical, it is not so weak that it has no dynami¬ 
cal effects. Indeed, our simulations give the counter-intuitive 
result that the Lorentz force in the initially weak-field simu¬ 
lation has a larger dynamical effect than that in the initially 
strong-field simulation. The density dependence of the mag¬ 
netic field in both molecular clumps and in our simulations 
appears to be an emergent property associated with turbu¬ 
lent, magnetized, self-gravitating systems. A more complete 
understanding must await further research. 

6.1.1 Scaling with Velocity Dispersion 

Mouschovias (1991) suggested that the magnetic field in 
molecular clumps depends on the velocity dispersion, cr„, as 
well as on the density, B oc Based on HI and OH Zee- 

man data from Crutcher (1999), a significantly smaller sam¬ 
ple than that in Crutcher et al. (2010), Basu (2000) showed 
that this magnetic field scaling is in excellent agreement with 
observation. He also found that the velocity dispersion is 
independent of density. Basu (2000) derived this relation 
under the assumption that clouds are planar and concluded 
that the agreement with observation supported this hypoth¬ 
esis. In Appendix A we show that this relation comes from 
an identity and applies to clouds of a variety of shapes pro¬ 
vided that they are gravitationally bound and are close to 
being magnetically critical, and provided that and avir 
are independent of density. Now with all 68 OH±CN sources 
from Crutcher et al. (2010), we confirm the prediction of 
Mouschovias (1991) and Basu (2000), Bioa/cry oc 
although there is a significant uncertainty in the exponent. 
However, as we have noted above, the magnetic field-density 
relation is in fact Bios oc implying that the veloc- 
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(i?tot/n?i®")/(Stot/n?i®')med (i?tot/n?i®")/(i?tot/n?i®'')med 

Figure 7. (a) The PDFs of normalized by the median value of the clump sample, from the strong-field model at 0.57tfi 

(shaded histogram) and from the OH-I-CN sources (open histogram) in Crutcher et al. (2010). K-S tests show a good match of the PDFs 
at this time and at different times as well, (b) Same as (a) for the weak-field model at 0.94tff. K-S tests show that the weak-field model 
agrees with observation about as well as the strong-field model, (c) The PDF of the normalized values of of the 100-clump 

sample from the strong-field model at 0.57tff compared with the uniform distribution proposed in Crutcher et al. (2010). (d) Same as 
(c) for the weak-field model at 0.94t(j. 


ity dispersion scales weakly with density. Indeed, a direct 
fit of (T„ vs. n gives CTi, oc as expected from these 

two relations. We also find that the observations imply that 
the valnes of /r* and Ovir inferred from line-of-sight data are 
independent of density to within the errors. 

The simulations give somewhat different results. For 
the strong-field model, we find that the time-averaged cr„ oc 
n®'35±0 03 ^ j.^o.i3±o.04. independent of density 

to within the errors. Despite the fact that the two power- 
law indexes increase slowly with time as the clumps undergo 
gravitational collapse, the value of a does not exhibit a sec¬ 
ular variation with time for t > O.Ats, although it does fluc¬ 
tuate somewhat (Fig. 6). The 100-clump sample in the weak 
field model behaves similarly, with small differences in the 
power-law indexes. The results of both models are consis¬ 
tent with Equation (A8). Observationally, we do not have 
information on the time evolution of these relations. Fit¬ 
ting these power-law relations to the OH and CN sources 
separately results in very large uncertainties because of the 
small number of sources and the limited range of densities 
in each group of sources. Therefore, although it is plausible 


that the CN sources, being denser, represent a later stage of 
evolution than the OH sources, current observations do not 
enable us to determine if the power-law indexes in actual 
molecular clumps evolve in time. 


6.2 Distribution of Magnetic Field Strengths 

We first show that distribution of line-of-sight fields in 
our simulations agrees remarkably well with observation. 
In Fig. 7a, we plot the probability density function (PDF) 
of normalized by the median (Hios/hH®3)med, of 

our 100-clump sample from the strong-held model at 0.57tfi 
together with that of the total 68 OH-I-CN sources from 
Crutcher et al. (2010). The PDF of the 100-clump sample 
is normalized to a total of 68 sources for direct compari¬ 
son. We performed a K-S test on these two PDFs, with the 
null hypothesis being that the two samples are from the 
same continuous distribution. The result is that the proba¬ 
bility that this could be true (the p-value of the test) is 0.89. 
We performed K-S tests on the simulation data at different 
times with different binning sizes. The tests gave p-values 
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Figure 8. Lognormal fit of the normalized values of of the 100-clump sample from (a) the strong-field model at 0.57tff and 

(b) the weak-field model at 0.94tff. K-S tests show a good match of the PDFs with p—values of 0.89 and 0.98, respectively. 


between 0.50 - 0.99 in the strong-held model. We have per¬ 
formed K-S tests on the weak-field model as well, and the 
results are similar to those for the strong field model. The 
comparison of the strong- and weak-field model results will 
be discussed in Section 8. We conclude that the PDFs of 

from the simulations and from 
the observed OH-I-CN sources are consistent with being 
drawn from the same continuous function. 

An important conclusion from the Bayesian study in 
Crutcher et al. (2010) is that the inferred distribution of 
the total field strength, Btot, is best described by a uniform 
distribution over the range f Btot, max < -Btot < Btot , max 
with / = 0.03 and Btot, max oc h)). Here we shall look at the 
shape of the distribution based on the values of Bios, Btot, 
and / from our simulations. We begin by considering how 
the observed and simulated line-of-sight fields agree with ei¬ 
ther a uniform or log-normal distribution. As we found in 
Section 3, a 'C test suggests that a uniform distribution for 
Btot does not provide a good fit to the observed line-of-sight 
data, since ~ 1.85. On the other hand, a K-S test, which 
is weaker, indicates that a uniform distribution is quite con¬ 
sistent with the observed data, since the p-value is 0.97. A 
log-normal distribution provides a slightly better fit to the 
observed line-of-sight data, with = 1.45; the p-value for 
the K-S test is 0.98. These results are inconclusive; the ob¬ 
served line-of-sight data does not allow us to infer whether 
the underlying distribution of total field strengths has a uni¬ 
form, log normal, or some other distribution. Analysis of the 
line-of-sight results from our two simulations yields a simi¬ 
lar result, although both a uniform distribution and a log¬ 
normal distribution provide better fits than they do to the 
observed data. 

However, the 3D data available from the simulation 
provides a much more powerful discriminant than the line- 


of-sight data: The distribution of the total field strength, 
Btot/h-H®®/(Btot/hH®®)med, from our simulations is not con¬ 
sistent with a uniform distribution-the null hypothesis is 
rejected with a p-value of 0.004 in the strong-field case and 
0.001 in the weak-field case. This is shown graphically in 
Fig. 7c and 7d. On the other hand, the distribution of nor¬ 
malized total field strengths from the simulation is well fit 
by a log-normal distribution (Figs. 8a,b): A K-S test gives 
a p-value of 0.89 for the strong-field case and 0.98 for the 
weak-field case, with o = 0.21, 0.28 for the two cases, re¬ 
spectively. The peak of the lognormal distribution is very 
close to the median value of the observed distribution in 
both cases. 

We conclude that both the observed and simulated line- 
of-sight data are consistent with each other and with a uni¬ 
form distribution of Btot. On the other hand, the simulated 
data for the total normalized field strength is not consis¬ 
tent with a uniform distribution, but it is consistent with a 
log-normal distribution. If we knew that at least one of our 
simulations was a reasonable representation of reality, we 
could conclude that observed fields do not have a uniform 
distribution of normalized total field strength. However, the 
many differences between the simulation and the observa¬ 
tions discussed in Section 5 prevent us from drawing such a 
conclusion. 

An important implication of the the work of Crutcher 
et al. (2010) is that magnetic field strengths in molecular 
clumps extend to very low values, Btot, min — 0.03Btot, max. 
From Equation (3), this implies min(Btot,/j/h-H^I) = 2.9. 
Now, the median value of Bios/uh®! = 11-9 from the data 
in Crutcher et al. (2010), so their results imply 

— ' ^ /-o ~ (Crutcher et al. 2010). (25) 
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In other words, they infer that the minimum normalized 
total field is about 1/4 of the observed median line-of-sight 
field. We can actually measure this ratio in our simulation, 
and in our strong-field case we find that this ratio is 0.69, 
almost 3 times as large. 

In addition, the total field is subject to cancellation, 
since it is an average over a Gaussian beam. The total field 
strength is better characterized by the rms field. Averaging 
over the whole clump (and not just the central beam as in 
the observations), we find 


min (I?rms, 4 )clump 

(Blos/h“®|)med 


= 1.07, 0.93 


(26) 


for the strong- and weak-field cases at 0.57 tg and 0.95 ta, re¬ 
spectively. (Note that the numerator is for the whole clump, 
whereas the denominator is the observed line-of-sight value 
for the central beam.) Thus we find no evidence for very 
weak fields in our simulation; instead, the minimum rms 
field is comparable to the median line-of-sight field. We also 
find that 


(Brms.v/Blos)med = 2.03, 1.94 (27) 

for the strong- and weak-field cases, respectively, where 
firms,V is the volume-weighted firms averaged over the vol¬ 
ume of the whole clump. Since fiios is mass-weighted and on 
average is half of fitot (also mass-weighted), it follows that 
firms.v is very nearly the same as the mass-weighted fitot 
averaged over a 0.025 pc radius central Gaussian beam. 


6.3 Possible Dependence of fiios/ng on Beam Size 

In Section 3, we pointed out that the 48 OH and GN sources 
for which we could determine the size of the telescope beam 
show a weak but significant dependence of fiios/h-H®® on 
rbeam. Fig. 3 shows that we see a similar effect in our 
strong-held simulation. The three solid circle symbols show 
(filos/h&®®)med for beam sizes of 0.025 pc (the one used in 
the analysis above), 0.0125 pc and 0.00625 pc. For our 2-level 
simulation, the results with a beam size of 0.00625 pc are 
not converged, so the value of (fiios/hH®®)med for this case is 
obtained from our 3-level run in the convergence study. The 
three data points from the simulations match almost exactly 
the best-htting slope of the 48 observed cloud clumps. Since 
hn almost certainly decreases with rseam, the increase of 
fiios /h-H®® with rbeam implies that the value of a is less than 
0.65 within individual clumps, possibly due to a diffusion 
mechanism operating within the clumps. Our simulations 
are based on ideal MHD, so the agreement between the sim¬ 
ulation results and observation suggests that, at least down 
to about 0.025 pc, non-ideal MHD diffusion mechanism such 
as ambipolar diffusion are not responsible for the reduction 
in fiios/h-H®® on small scales. A possible candidate for this 
is reconnection in a turbulent medium (sometimes termed 
‘reconnection diffusion’—see Lazarian et al. (2015) and ref¬ 
erences therein), which will be discussed in Section 7. How¬ 
ever, because of the large scatter of the observed data in Fig. 
3, we cannot draw a firm conclusion on this possible beam- 
size effect. Future high resolution Zeeman observations will 
provide more information on this issue. 


7 OBSERVED VS. ACTUAL MASS-TO-FLUX 
RATIO 

As discussed in Appendix A, the mass-to-flux ratio, /r<t, is 
of fundamental importance in star formation since it deter¬ 
mines whether gravity can overcome magnetic stresses and 
initiate gravitational collapse. The observed mass-to-flux ra¬ 
tio is based on the mass-weighted line-of-sight field, fiios, 
and the column density, Ebeam, in the central beam: 



For the moment, let us assume that the clump is static. 
There are then two main corrections that are needed in or¬ 
der to obtain the actual mass-to-flux ratio: First, the field 
and column density must be measured along the direction of 
the mean field, which we label z'. This de-projected value, 
termed the ‘intrinsic value’ by Heiles & Grutcher (2005), is 
(fitot) in our notation. Second, the magnetic field must 
be changed from fitot, the mass-weighted field in a central 
Gaussian beam, to fimid.z', the field in, and normal to, the 
midplane of the clump with no mass weighting. The method 
for determining the midplane in a non-uniform clump will 
be described in Section 7.2 below. For a static, equilibrium 
cloud embedded in a uniform magnetic field, fimid, z' is the 
total field. We define fimid. z' so that it is averaged over a 
disc of radius rbeam at the center of the clump; this is prefer¬ 
able to averaging over the whole clump since the size of the 
clump depends on the somewhat arbitrary criterion adopted 
to define the clump. In addition to correcting the field, the 
column density must also be corrected so that it refers to 
the mass inside a cylinder of radius rbeam rather than that in 
a Gaussian beam of that characteristic radius; we label this 
Elmid.z'- The result, /r* (fimid, z'): i® the actual mass-to-flux 
ratio for a static clump. 

In the real world, clumps are not static, and the field 
can become tangled. The critical value of the mass-to-flux 
ratio is not known in this case; we shall continue to use the 
value for a static cloud here. 


7.1 Projection Effects 


The distribution of the observed field, which is the line-of- 
sight field, to the distribution of the intrinsic field is given 
by Equation (B5). One can immediately show that this im¬ 
plies that (1/fiios) is infinite. As a result, it is preferable 
to work with the flux-to-mass ratio, ^ fiios, in eval¬ 

uating the average properties of the observed mass-to-flux 
ratio. Even when dealing with the intrinsic mass-to-flux ra¬ 
tio, which does not suffer from this divergence, it is better 
to average if the total field can assume very small val¬ 
ues, as inferred by Crutcher et al. (2010). The mean of the 
flux-to-mass ratio is the inverse of the harmonic mean of the 
mass-to-flux ratio. 


(flux-to-mass ratio) 



1 


(29) 


For spherical clouds, the column density is not affected 
by projection, so that (/i*,ios)^ = 2(/r,t)^ since (fiios) = 
5 (fitot). The median values of fiios and /i<t,ios depend on the 
distribution of field strengths, 4>{Btot)', for a delta-function 
distribution, fiios,med = |fitot and ^.j,ios,med = 2^*(fitot). 
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Table 4. Line-of-Sight and Actual Mass-to-Flux Ratios of the 100-Clump Samples 


Model“ 

t/iff 

II ^ 

A^<i>,los 


//^(Smid, 2' 

')° 

/^«(|B|mid, 

z>r 



H. Mean 

Med 

H. Mean 

Med 

H. Mean 

Med 

Strong 

0 

1.32 ±0.10 

1.65 

0.95 ± .05 

0.86 

0.91 ± .05 

0.81 

field 

0.4 

2.00 ±0.17 

2.53 

1.57 ± 0.09 

1.69 

1.52 ±0.08 

1.67 


0.5 

1.95 ±0.19 

2.63 

1.56 ±0.10 

1.83 

1.51 ± 0.10 

1.70 


0.57 

2.27 ± 0.18 

2.87 

1.80 ±0.12 

2.05 

1.68 ±0.10 

1.94 


0.64 

2.25 ±0.19 

2.96 

1.75 ±0.13 

1.98 

1.63 ±0.11 

1.86 


Time avg.® 

2.08 ±0.10 

2.72 

1.63 ±0.06 

1.86 

1.56 ±0.05 

1.77 

Weak 

0 

1.71 ± 0.14 

2.36 

0.99 ± 0.05 

1.02 

0.92 ± 0.04 

0.97 

field 

0.5 

2.01 ± 0.18 

2.62 

1.32 ±0.12 

1.48 

1.25 ±0.11 

1.43 


0.6 

2.15 ±0.23 

3.19 

1.55 ±0.11 

1.82 

1.45 ±0.10 

1.72 


0.7 

2.22 ±0.23 

3.26 

1.52 ±0.14 

1.92 

1.42 ±0.12 

1.70 


0.8 

2.61 ± 0.27 

3.63 

1.99 ±0.23 

3.09 

1.82 ±0.19 

2.70 


0.9 

2.66 ±0.24 

3.51 

1.83 ±0.20 

2.76 

1.66 ±0.19 

2.37 


0.94 

2.57 ± 0.29 

4.00 

2.04 ±0.20 

3.02 

1.88 ±0.18 

2.51 


Time avg. 

2.29 ±0.10 

3.30 

1.56 ±0.06 

2.27 

1.45 ±0.06 

2.02 

Observed-^ 


2.40 ±0.29 

4.10 

- 

- 

- 

- 


^ See Table 2. 

Observed ^ 13.4051 which is derived from the line-of-sight column density and -Bios by convolving the clump with a Gaussian beam of 
0.025 pc in radius (see Section 6). 

Actual which is based on the field normal to the clump midplane (see Section 7). 

Value of /i<j based on the magnitude of the field normal to the clump midplane (see Section 7). 

® The time average excludes t = 0. 

^ From the 68 OH-I-CN clumps in Crutcher et al. (2010) 


At the other extreme, if clouds are flattened into sheets, 
then the observed column density is larger than that nor¬ 
mal to the sheet by a factor 1 /cos 9, so that the average of 
the line-of-sight values is related to that of the total-field 
values by {n<s>,ioB)h ~ (this is a corrected version 

of the result given in Heiles & Crutcher 2005); for a delta- 
function distribution of held strengths, the median is given 
by /i<s,ios,med = d/T,} (.Btot)• Li et al. (2011) have pointed out 
that gravitational contraction along field lines can increase 
the apparent mass-to-flux ratio, but the maximum possible 
value In their model Is that for a sheet, so their discussion 
fits within the existing framework. The medians and har¬ 
monic means of /i$,ios of the 100-clump samples at different 
times from the strong and weak held models are given in 
Table 4. 


7.2 The Actual Mass-to-flux Ratio 

The clumps in a turbulent medinm do not have the axisym- 
metric shape assumed in theoretical treatments of the mass- 
to-flux ratio. Lacking an axis of symmetry, along what di¬ 
rection do we evaluate /rs? Imagine evaluating along dif¬ 
ferent directions in the axisymmetric case. Along the axis of 
symmetry, the flux will be a maximum and the column den¬ 
sity a minimum, so that will be a minimum. Requiring 
the minimum to exceed unity provides the most strin¬ 
gent test for whether the clump is supercritical. We therefore 
adopt this criterion in the turbulent case. To implement it, 
we divide the polar angle into 24 divisions (7.5° apart) and 
the azimuthal angle into 48 divisions and evaluate in the 
central at each orientation. The orientation of the 

midplane is the one that gives the minimum value of 


the center of the midplane is defined as the cell with the 
maximum density, and the midplane itself Is defined by the 
condition that it pass through this cell. 

To obtain the actual mass-to-flux ratio, we take the 
steps described above. The ratio of the actual mass-to-flux 
ratio to the line-of-sight value is 

actual _ Ljniid, z' Rios /ont 

II T ' R ' Ibb) 

los -^beam -^mid, z' 

Recall that .Bmid, z' is the average value of the component of 
the field normal to the midplane over an area TTr^eam s-t the 
center of the clump; it is not density weighted. The line-of- 
sight field, -Bios, differs from this in several respects: It is den¬ 
sity weighted, it is averaged along the line of sight, and it is 
also averaged over a Gaussian beam. One typically assumes 
that the line-of-sight field is about half the actual field due 
to projection effects. In fact, the density weighting of the 
line-of-sight field increases the ratio; we find that the me¬ 
dian value of Bios/Bmid 2 ' is 0.80. As noted above, spherical 
clouds have isotropic column densities, whereas sheets lead 
to an increase in the ratio of the observed column density to 
the value that enters For the clumps in our simulation, 
the median value of the ratio of the actual column density 
to the line-of-sight value (which is averaged over a Gaussian 
beam) is 1.09, close to the spherical value. Equation (30) 
then implies 

actual, med — 0. 7/r4>, los, med . (31) 

The key result of this analysis is that the actual mass-to- 
flux ratio is significantly greater than (0.25-0.5) times the 
observable line-of-sight value, which is what is expected from 
a naive application of the projection effects discussed above. 
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Detailed results for the mass-to-flux ratio are given in 
Table 4. The ratio of the median of the actual /r* to the 
median of the line-of-sight value is 0 . 68 , very close to the 
estimate above. The harmonic means are closer: the actual 
value is 0.78 times the line-of-sight value. The weak mean- 
field results are similar to the strong mean-field ones. 

It should be noted that clump mergers and the growth 
of clumps due to flows along field lines imply that the iden¬ 
tity of the 100 clumps in the samples can evolve with time. 
Mergers of clumps in the samples allow new, smaller clumps 
with lower mass-to-flux ratios to enter the samples; field- 
aligned accretion leads to an increase in the mass-to-flux 
ratio. The latter effect dominates, so that there is a slow 
increase in the mean mass-to-flux ratio with time. 

Next, we assess the effects of tangling of the field by 
evaluating /r* with the magnitude of the held in the mid¬ 
plane, |B|mid, z', rather than the component normal to the 
midplane, Bmid, z'; the results are only slightly smaller, 
showing that held-line reversals are insignihcant in the mid¬ 
plane (the time-averaged median value of |i3mid, z'|/^mid, z' 
is about 1.01). Field reversals are somewhat more impor¬ 
tant for the line-of-sight held, which averages over a volume 
along the line of sight: The time-averaged median value of 
(|Bios|/Bios) for the 100 clumps is 1.29. As shown in Ta¬ 
ble 4, comparison of the time-averaged /r# (Bmid, z') with 
z') shows that the ehect of held-line tangling in 
the weak-held case is a little larger than that in the strong- 
held case, which is a natural outcome for the weak held 
model. 

The harmonic mean of p<j,ios from the observed 68 
OH-bCN clumps in Crutcher et al. (2010) is about 2.4 as 
shown in Table 4, close to the values for our 100-clump 
samples. (We do not place great weight on this agreement, 
however, since we have found that the values of in our 
samples depend weakly on mass, an ehect we shall explore 
in a subsequent paper.) The results of our strong-held sim¬ 
ulation suggest that the harmonic mean of the actual mass 
to hux ratio is about 78 per cent of the line-of-sight value, 
or actual,/i — 1.5. However, because the harmonic mean 
is disproportionately ahected by small values, the median 
value is more characteristic of the observed mass-to-hux ra¬ 
tio. Since the observed median is /i<j,ios.med — 4, we infer 
that the actual median mass-to-hux ratio is /i 0 , actual, med — 
2.8 from Equation (31). Typical clumps in molecular clouds 
are therefore magnetically supercritical by a signihcant mar¬ 
gin (~ factor 3). 

A striking result of the strong-held simulation is that 
the time-averaged median of the mass-to-hux ratio in the 
100 most massive clumps, WiiBmid, z') = 1 . 86 , is somewhat 
greater than the initial value for the entire simulation box, 
P<t.o = 1.62. In ideal MHD, the only way in which the 
mass-to-hux ratio can increase beyond its initial value is 
if there are held-line reversals in the plane used to calculate 
the hux. However, although the time-averaged median of 
p<t(|i3|j„id, z')) which excludes the held-line reversal effect in¬ 
side the center of the clumps, is smaller than time-averaged 
median of /^^.(Bmid,z'): il larger than /i 4 .,o. Fully 25 

per cent of the clumps have mass-to-hux ratios based on 
|S|mid,z' that are at least twice the initial value. This can 
occur only due to numerical ehects. These are of two types: 
One is purely numerical and can be mitigated by increasing 
the resolution. The convergence study in the Appendix indi¬ 


cates that the strong-held simulation is reasonably well con¬ 
verged, which implies that this type of numerical ehect is not 
important. The second type is the numerical representation 
of microphysical processes such as viscosity and reconnec¬ 
tion. For example, simulations of turbulence often rely on 
numerical viscosity to dissipate energy on small scales; as 
the resolution is increased, the scale at which the energy is 
dissipated shrinks, but (provided there is a well-established 
inertial range) the rate of dissipation does not change. Since 
the increase in the mass-to-hux ratio above the initial value 
is converged, the numerical ehect that allows this is of the 
second type. A process that is known to violate hux-freezing 
is turbulent reconnection (see the review by Lazarian et al. 
2015), and our results on the increase in /r*(|B|mid, z') over 
its initial value are consistent with a numerically enabled 
version of this process. 

Finally, we note that Bertram et al. (2012) have ana¬ 
lyzed the spatial distribution of 77. = /i*, core/M*, envelope in 
simulations of turbulent cores. They found that this ratio 
is generally > 1 , although there is a large dispersion with a 
non-negligible fraction of cases with 77 < 1. A value 77 > 1 is 
expected in theories of star formation that rely on ambipo- 
lar dihusion Crutcher et al. (2009), but the simulations of 
Bertram et al. (2012) were based on ideal MHD. (Crutcher 
et al. 2009) found 77 < 1 in each of the four sources they 
observed; Bertram et al. (2012) attribute the difference be¬ 
tween their simulations and the observed results to the small 
size of the observed sample. Our simulations have a higher 
resolution than those of Bertram et al. (2012) (512^ base 
grid plus two levels of refinement vs. 256®) and had a longer 
driving period (2 crossing times vs. 1 crossing time), but 
we find a similar result: the value of for entire clumps 
(cores plus envelopes) is, on average, greater than that for 
the central beam. Since our results are consistent with those 
of Bertram et al. (2012), we have not carried out a detailed 
analysis of the spatial dependence of the mass-to-flux ratio. 

7.3 Comparison with CN Zeeman Observations 

Whereas the OH observations come from gas with densities 
less than that of the clumps in our 100 -clump sample, the 
CN observations (Falgarone et al. 2008) are from gas with 
densities that lie within the simulated range, so we consider 
these sources separately. In fact, only 16 of our 100-clump 
sample have both a density and a column density that is in 
the observed range for CN clumps (hn > 2x 10® cm“®, TVh > 
10 ®® cm“®), so we create subsamples with these properties. 
In Fig. 9, we plot the line-of-sight quantities Bios and /i 4 .,ios 
versus Nu for these subsamples together with the results 
for the 14 observed CN sources. Since the median beam 
size for the CN sources is comparable to the size of the 
sources, we compare the observed values of Bios and /r$,ios 
with simulated values from the whole clump instead of the 
smaller central-beam region. We can see that the simulated 
clumps are basically in the same region occupied by the CN 
sources, with a few more clumps from the weak field model 
lying on the smaller Bios part of Fig. 9a and larger /i*,ios part 
of Fig. 9b compared to the strong field model. Again, both 
models reproduce the location of the observed CN sources 
in Fig. 9. There is one CN source, out of the 14 observed, 
that is sub-critical based on the observed value of Bios, and 
two others that are potentially sub-critical since Bios is on 
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Figure 9. Comparison of results from the strong-field simulation with observations of CN clumps (Falgarone et al. 2008), which have 
A^h > 1 X 10^® cm“^ and nn > 2 X 10® cm“®. The 14 observed CN clumps are labeled with green solid circles; 1/3 of the lines of 
sight from the strong-field 100-clump sample at 0.57iff are labeled with black circles; and 1/3 of the lines of sight from the weak-field 
100-clump sample at 0.944ff are labeled with blue squares. Since the beam size for the CN observations is comparable to the clump size, 
we have calculated the field strength and column density by averaging over the entire clump for this figure only. The distribution of (a) 
Bias S'lid (t>) M*,los of clumps from both models matches well with data for the CN sources. The straight lines in the panels indicate the 
magnetic critical condition. See Section 7.3 for discussion. 



Figure 10. Streamline plots of the magnetic field in the two simulations at t = 0. The difference is clear between the relatively ordered 
magnetic field in the strong-field model (a) and the highly perturbed and stretched field in the weak-field model (b). The highly perturbed 
field in the weak-field model makes it more difficult for material to flow along along the field lines. 


or near the critical lines in Fig. 9. One of these two CN 
sources, S106OH, would be subcritical if the line-of-sight 
field is less than 1 sigma below the measured value. The 
other one, G10.6, would be subcritical if the actual mean 
field direction makes an angle larger than 24° from the line 
of sight. Interestingly, it is the weak-field model that has 


a sub-critical clump and more clumps close to the critical 
line. In general, the distribution of Bios versus A^h of the 
high density clumps from our simulation matches well with 
that of the CN sources. Future observation using ALMA will 
be able to increase the number of CN sources for a better 
statistical comparison with our simulation. 
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8 IS THE MEAN FIELD STRONG OR WEAK? 

8.1 Comparison of the Strong- and Weak-Field 
Models 

The global properties of the strong- and weak-field simula¬ 
tions are shown in Table 2. Initially the uniform magnetic 
fields differ in strength by an order of magnitude. After tur¬ 
bulent driving for two crossing times, the magnetic held in 
the weak-held model is enhanced enough to increase firms,v 
by more than a factor 5, so that the Alfven Mach num¬ 
ber drops from Ma ~ 10 to A4a = 1.86. By contrast, in the 
strong-held simulation. At a drops by only about 15 per cent. 
These results are portrayed graphically in Fig. 10, which 
shows the held lines of the magnetic held inside the turbu¬ 
lent boxes. The magnetic held is roughly along the mean 
held direction in the strong-held model (Fig. 10a), but it is 
highly tangled and stretched in the weak-held model (Fig. 
10b). Initially, the values of firms, 4 of Ih® systems 
diher by a factor of 10, but after two crossing times, the 
diherence is less than a factor of 2. At the end of the simula¬ 
tions, the rms magnetic held has been enhanced by a few per 
cent (strong held) or slightly more than 10 per cent (weak 
held) due to gravitational contraction. 

Next, consider the results for the 100-clump samples. As 
shown in Table 3 the magnetic-held properties of the clumps 
in the strong-held model provide a good match with ob¬ 
servation. One might therefore expect the weak-held model 
to have results that clearly differ from observation, since it 
has an initial held strength that is a tenth of that in the 
strong-held model and, correspondingly, a plasma /3 that is 
a hundred times that of the strong-held model (Table 2). 
We identify the 100-clump sample in the weak-held model 
in the same manner as in the strong-held model, with the 
same density threshold of hu = 4 x 10^ cm“®. In Fig. 11, we 
repeat Fig. 1 and then include the data points from the 100- 
clump weak-held sample at 0.94tff. The data points from the 
two models overlap substantially and are not readily distin¬ 
guished from each other. 

The time-averaged value of the power-law index a(fitot) 
between fitot and hu for the weak-held model is 0.57T0.05, 
signihcantly less than the value inferred from observation, 
O.GS^Q jg. The strong-held model has a signihcantly larger 
value, a(fitot) = 0.70T0.06, but it agrees with the observed 
value within the errors. One would have expected the weak- 
held model to have a density dependence closer to the value 
of 2/3 predicted by Mestel (1966) for weak helds. However, 
as we have noted above, it is only the mean held that is weak 
in the weak-held model; turbulence has amplihed the small 
scale held so that it is comparable to that in the strong-held 
model. Furthermore, the tangling of the held can be ampli¬ 
hed as the clumps contract under the inhuence of gravity. 
We have conhrmed this picture by calculating (Q(|fiios|))j 
for the magnitude of the line-of-sight held. For the strong- 
held case, this is 0.67 ± 0.04, which is consistent with the 
value for the observable line-of-sight held given in Table 3, 
(a(Bios))t = 0.63 ± 0.07. By contrast, for the weak-held 
case, we hnd that (a(|fiios|))^ = 0.64 ± 0.03 for the magni¬ 
tude of the line-of-sight held, whereas Table 3 shows that 
(a(fiios))t = 0.58 ± 0.05. We conclude that the magnitude 
of the held increases with density at about the same rate for 
both the strong and weak-held cases (consistent with the 



Figure 11. Results of the weak-field model at 0.94 tg (blue 
crosses) are overplotted on the data in Fig. 1. Both models pro¬ 
duce clumps in the same region of the Bios vs nn diagram, but 
the median value of Bios in the weak-field model is 80 per cent of 
that in the strong-field model. 


observed value, a(fiios) = 0.64 ± 0.13), but that the rate 
of increase for the observable line-of-sight held, fiios, is less, 
just as it is for fitot. The increased tangling of the held with 
compression in the weak-held case prevents the observable 
held from increasing at the expected rate of a = 2/3. 

Table 3 also gives the means and medians of 
fitot, ti/h-H®! and fiios, ti/h-H,®! for the weak-held sample. 
These quantities are plotted in Fig. 6d-e as functions of time 
during the gravitational collapse phase; they are relatively 
constant for t > 0.5ta. The time-averaged means and medi¬ 
ans of these quantities agree with the strong-held values to 
within about 25 per cent, even though the initial helds diher 
by a factor 10. The mean and median values of p#,ios for the 
weak-held model are plotted as a function of time in Fig. 6f; 
they are larger than those in the strong-held model by about 
a factor of 1.4, and after t = 0.5tg they increase only slightly 
over the remainder of the simulation. In Fig. 12, we plot hn 
and fitot of the 100-clump samples from the two models at 
t — 0. The mean magnetic held in the weak-held clumps, 
(fitot)t, is about 70 per cent of (fitot)t of the strong-held 
model clumps, and the median densities are even closer. 

Finally, we compare the results for the 100-clump sam¬ 
ples with the results for the entire system; in this case the 
held in the clumps is volume-averaged rather than mass- 
weighted, since it is the volume-averaged held that is rele¬ 
vant for clump structure. As we discussed above, the rms 
held changes in the turbulent driving phase much more for 
the weak-held case than the strong-held one. Including the 
subsequent phase in which self-gravity is operative, we see 
from Table 2 that from the initial uniform held to the end 
of the simulations, firms, v of the whole system increases by 
a factor of 1.2 in the strong-held case and about a factor 
of 6 in the weak-held model, mainly due to enhancement of 
the held by supersonic turbulence. On the other hand, the 
median values of firms, v in clumps in the 100-clump sam¬ 
ple are about 2.7 and 4.6 times the values of firms, v for the 
entire turbulent box for the strong- and weak-held models. 
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Figure 12. Distribution of Btot vs njj in a 0.025 pc radius Gaus¬ 
sian beam for the 100-clump samples from the strong-held (black 
circles) and the weak-Held (blue squares) models at t = 0. The 
median values of Btot and nn of the strong-held sample (solid 
lines) and the weak-held sample (dashed lines) are also shown. 
The clumps from the weak-held model occupy the low-Btot re¬ 
gion of the plot. 



Figure 13. The mass fraction, /jv/, of gas with density njj ^ 
1 X 10^ cm“® (thin blue lines) and riH ^ 4 X 10^ cm“® (thick 
red lines) of the strong (solid line) and weak (dashed line) mean 
held simulations as a function of time. Material can be accreted 
onto dense clumps more easily with the ordered magnetic Held in 
the strong-held model than with the highly perturbed held in the 
weak-held model (see Fig. 10). The difference in accretion rates 
onto the dense regions of the two models is about a factor of 1.5. 


respectively, at the end of the simulations. The higher val¬ 
ues of firms,V in the clumps compared to that in the entire 
simulation box is simply the result of further enhancement 
of the held by the gravitational collapse of dense clumps. 

The remarkable similarity in these results of the strong- 
and weak-held simulations is due to the greater enhancement 


of the magnetic held in the weak-held model as the result of 
both easier compression and greater stretching of the held 
lines. Correspondingly, the turbulent Alfven Mach number 
in the weak-held model has dropped from 10 to 1.61, not 
much larger than the the hnal value in the strong-held case, 
0.83, at the end of the simulations. Similar results have been 
found previously (e.g. Lunttila et al. 2008; Padoan et al. 
2010). However, there are two important diherences between 
the two models. Prior to the time that gravity is turned 
on, the median mass-to-hux ratio does not change signih- 
cantly in the strong-held case, but it drops by more than 
a factor 15 in the weak-held case. This is consistent with 
turbulent fragmentation along hux tubes: the median held 
in the weak-held simulation has increased by a factor 6.7 
at t = 0, and the median surface density in the 100-clump 
sample is 0.58So. As a result, oc E/fi should typically 
drop by a factor of about 11.6. This is consistent with the 
results in Table 4, which gives z')med = 15.9, 

about 1/3 larger than our simple estimate. Note that by us¬ 
ing the mean held in the box rather than in the clumps in 
our estimate, we implicitly assumed that the clumps formed 
primarily by how along held lines rather than by compres¬ 
sion. The agreement with the simulation results shows that 
this is approximately correct. 

We can develop a simple picture of the turbulent frag¬ 
mentation model as follows. In the strong-held case, the fact 
that the mass-to-hux ratio of the most massive clumps (i.e., 
the ones in our 100-clump sample) is typically about the 
same as the initial value implies that these clumps form by 
accreting material along a hux tube over a distance compa¬ 
rable to the size of the box. What about the weak-held case? 
In this case, the held is stretched and compressed. So long as 
hux-freezing holds, we can dehne hux tubes with dimensions 
that evolve smoothly in time. Since the volume in the sim¬ 
ulation is hxed, the average volume of a hux tube must be 
constant. Consider a hux tube in the initial uniform held, 
with a length £o equal to the box size, an area normal to 
the held Aft,o, and a volume Vft,o = A{tfitit,o- After turbu¬ 
lence begins, the average area of a hux tube shrinks by the 
same factor as the length increases. Aft = Vtt.o/fft- Conser¬ 
vation of magnetic hux then implies fi = fio(Aft,o/Aft) = 
Bo{£ft/£o)- the held increases in direct proportion to the 
length of the hux tube. Mass conservation gives a similar 
result, E = Eo(£ft/f'o), so that the mass-to-hux ratio E/fi is 
constant in time. Fragmentation along the hux tube reduces 
the mass and surface density by a factor ^frag/^ft, where 
£frag is the length from which the fragment mass is drawn, 
but fragmentation does not ahect the held. As a result, the 
mass-to-hux ratio in the fragment is (^frag/^ft)M 4 >,o. For ex¬ 
ample, the results of the strong-held simulation discussed 
above show that the typical hux tube was stretched by a fac¬ 
tor 6.7, increasing fi and E by the same factor, while leaving 
the mass-to-hux unchanged. Typical fragments were drawn 
from a length ^uag = 0.58^o- The typical mass-to-hux ratio 
was then reduced by a factor ^nag/^ft = 0.58/6.7 = 1/11.6. 

Once gravity is turned on, a second important diher- 
ence between the strong- and weak-held cases emerges: Ma¬ 
terial accreting onto the clumps is slower in the weak-held 
case than in the strong-held one because of strong magnetic 
held tangling. Fig. 13 shows that the rate of gravitational 
collapse, as measured by the rate of increase of the mass 
fractions above densities nn ^ 10'* and is 4 x 10* cm“^. 
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is about 1.5 times greater in the strong-field case than in 
the weak-field one. (Note that the factor ~ 1.5 difference 
between the rate of increase of dense gas in the strong-field 
model compared to the weak-field one does not mean that 
the star formation efficiencies must differ by a factor 1.5 
since the star formation rate is affected by feedback pro¬ 
cesses that we have not considered in this work.) We inter¬ 
pret this counter-intuitive result as being due to the much 
greater twisting and stretching of the field in the weak-field 
case (Fig. 10), which makes it more difficult for material 
to flow along the field and condense onto dense regions. At 
t = 0, before gravity has any effect, the mass fraction of 
dense gas, /m, in the strong-field model is already higher 
than that in the weak-field model. Since the turbulent driv¬ 
ing is the same in the two models, the only differences in 
the models are the strength and tangling of the field. In 
the strong-field model, compression along relatively straight 
field lines can create more dense regions than in the weak- 
field model, in which the more isotropic magnetic forces from 
the highly tangled field better resist the formation of dense 
regions. During gravitational collapse, the rate of increase of 
/m increases slightly in the strong-field model but remains 
about the same or even decreases slightly in the weak-field 
model. The median values of the ratio of the infall velocity 
to the free-fall velocity (Equation 23) of the 100 clumps in 
the weak-field model at 0.6 and O.Odtff are 0.17 and 0.15, 
respectively, about half the value in the strong-field model 
(0.29), a direct indication of slower material accretion onto 
the clumps. We compute the degree to which the Lorentz 
force, Fl, in each clump resists self gravity by summing up 
the component of the force opposite to the direction of self- 
gravitational acceleration, —Fl • g^ cell-by-cell, where g' is 
the unit vector associated with the gravitational acceleration 
due to the clump. Normalizing this to the total gravitational 
force, we evaluate 


LumpPS'-g'dV 


(32) 


for each clump. We find that the mean of this quantity for 
the 100-clump sample of the weak-field model is 1.35 times 
that of the strong-field model. This is a quantitative mani¬ 
festation of how the greater tangling of the field in clumps 
can inhibit gravitational collapse in the weak-field model. 

Table 4 lists the harmonic means and medians of /i* 
for the two models. At t = 0, after the initial turbulent 
driving, the median of the weak-field model is about 1.2 
times the value of the strong-field model, instead of the fac¬ 
tor of 10 from the initial conditions. At that time, more 
than half the clumps in the strong-field sample are sub- 
critical, whereas the majority of the clumps in the weak-field 
model are marginally supercritical. These results demon¬ 
strate that turbulent fragmentation is effective in reducing 
/i<s, in clumps, particularly in the weak-field case. 


8.2 Magnetic Field Orientation 

Recent polarimetric investigations of the orientation of mag¬ 
netic fields in cloud clumps on scales > 0.04 pc (Li et al. 
2009, 2011) show that it is not random, but instead is aligned 
with the global magnetic field of the surrounding intercloud 
medium (ICM, the medium between molecular clouds) on 


Table 5. Large-Scale Orientation of Magnetic Field 



% < 26° 

% < 45° 

Observed 

68 

84 

Strong-field simulation 
100 clumps 

25 clumps 

54 (56)“ 

70 (67) 

81 (80) 

86 (91) 

Weak-field simulation 

100 clumps 

25 clumps 

24 (24) 

22 (35) 

45 (42) 

42 (51) 

Random 

29 

50 


“ Numbers inside brackets are the proportions of magnetic field 
orientation obtained by integrating along the line of sight across 
the whole turbulent box assuming that it is optically thin. See 
Section 8.2 for discussion. 


scales of 100-200 pc. Li et al. (2009) found that almost 70 
per cent of the 25 clumps they observed have magnetic fields 
oriented within 26° of the ICM field, and that 84 per cent 
have fields oriented within 45° of the ICM field. As shown 
in Table 5, these percentages are much larger than expected 
by chance, and Li et al. (2009) estimate that the proba¬ 
bility that the 45° alignment occurs by chance is less than 
7 X 10“®. Comparing with several ideal MHD simulations, 
they concluded that only sub-Alfvenic models-i.e., models in 
which the magnetic field is dynamically significant in molec¬ 
ular clouds-can have such a significant correlation between 
the field orientation in molecular clumps and that in the 
ambient ICM. 

In Fig. 14, we plot the distributions of the magnetic 
field orientation with respect to the global mean field of the 
turbulent box for the two simulations. The global mean field 
is equivalent to the ICM field since it represents the field ori¬ 
entation of the material out of which the cloud was formed. 
Because the direction of the field in the plane of the sky is 
unknown, the magnitude of the relative orientation is in the 
range 0° — 90°. We plot the 200 relative orientations based 
on observations of the 100 clumps along the x- and y-axes, 
normal to the global mean field. In Fig. 14a, we plot the an¬ 
gle differences of the 100-clump samples of the strong- and 
weak-field models. The histogram is normalized to unity for 
direct comparison. As shown in Table 5, the strong-field sim¬ 
ulation shows alignments much greater than expected from 
chance, whereas the weak-field simulation shows no align¬ 
ment. In fact, a KS test shows that the weak-field results 
are consistent with a random distribution of orientations 
with p = 0.31. We can change the viewing direction so that 
it is not along the axes normal to the mean field, but this 
changes the percentage of clumps with alignments less than 
26° and 45° by less than 3 per cent. 

Most of the clumps in the 25-clump sample in Li et al. 
(2009) have sizes ^ 0.3 pc, much larger than the clumps in 
our simulations. If we focus on the largest 25 clumps in the 
strong-field simulation (Fig. 14b), which have sizes of 0.1- 
0.13 pc, we find much better agreement with the results of 
Li et al. (2009). Higher resolution observations in the future 
will be able to probe the magnetic field orientation closer 
to the center of the cloud clumps and determine whether 
fields on those scales are less aligned with the ICM fields. 
An initial step in this direction has been taken by Hull et 
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angle difference (°) angle difference (°) 

Figure 14. Orientation of the projected 3D, mass-weighted magnetic field, Btoti relative to the global mean field in the simulations, 
(a) The PDF of field orientations from the 100-clump sample in the strong-field (shaded histogram) and weak-field (open histogram) 
simulations. For the strong-field simulation, 54 per cent of the clumps have orientations within 26° of the mean and 81 per cent have 
orientations within 45°. The PDF for the weak-field simulation is consistent with a uniform distribution, (b) Same as (a) but only for 
the largest 25 clumps with nj] > Ucrit- Foi' the strong-field simulation, 70 per cent of the clumps have fields oriented within 26° of the 
mean, and 86 per cent are oriented within 45°, close to the observational results of Li et al. (2009). See discussion in Section 8.2. 


al. (2013), who find alignment within 45° between angular 
scales of 2.5" and 20" of 70 per cent. This is less than that 
found by Li et al. (2009), but the errors in the measurement 
of the orientations are large enough that this difference is 
not statistically significant. 

The above analysis includes only the material that we 
identified as being within a clump. Observers do not have 
this luxury; they must integrate along the entire line of sight. 
As an approximation to what observers do, we have obtained 
the mass-weighted field orientation by integrating through 
the turbulent box over the projected area of each clump 
in both the x and y directions. Clumps with mean column 
densities less than half the total are not included. This elim¬ 
inates 76 of the 200 clump observations in the strong-field 
model. There are two other differences between our method 
and that of observers: We did not attempt to use velocity 
information to separate out the gas that is associated with 
an individual clump; and we integrated only through the 
turbulent box, not along the entire line of sight through the 
clump. The results for the proportions of clumps with mag¬ 
netic field orientation within 26° and 45° of the mean-field 
direction are listed in Table 5 inside the parentheses. For the 
strong-field model, the new proportions are similar to those 
when the orientation of the magnetic field is computed only 
over the clumps; for the 25-clump sample, the proportion of 
sources varies by 3-5 per cent. For the weak-field simulation, 
the alignment increases somewhat in the 25-clump sample, 
but in all cases remains consistent with a random orienta¬ 
tion. Based on these results, one would expect that the line- 
of-sight field along the mean-field direction is larger than 
that normal to that direction for the strong-field case but 
not for the weak-field case, and this is indeed the case: The 
time-averaged median value of Bios,||/Bios,x = (1.86, 0.98) 
for the two cases, respectively; the weak-field case is consis¬ 
tent with isotropy. 

We conclude that the weak-field simulation is inconsis¬ 


tent with the observations of the alignment of fields on large 
and small scales by Li et al. (2009). The strong-field simula¬ 
tion has A4a ~ 1 and exhibits slightly less alignment than 
observed, consistent with the motions in observed molecular 
clouds being slightly sub-Alfvenic , as suggested by Li et al. 
(2009). It is also possible that in the strong-field case, am- 
bipolar diffusion could bring the small-scale field into closer 
alignment with the large-scale one, allowing a moderate in¬ 
crease in AIa. However, ambipolar diffusion would be in¬ 
effective in aligning the field in the weak-field case since 
ambipolar diffusion is ineffective in weak fields (e.g., McKee 
et al. 2010). 


9 CONCLUSIONS 

The measurement of interstellar magnetic fields is very chal¬ 
lenging: In the data compiled by Crutcher et al. (2010), 
there were fewer than 140 interstellar clouds with line-of- 
sight fields determined by Zeeman splitting, and among 
the molecular sources (about half the sample), only about 
40 per cent had fields detected with a significance greater 
than 2 sigma. Crutcher et al. (2010) carried out a sophisti¬ 
cated Bayesian analysis of this data, including the sources 
from which only upper limits are available, in order to in¬ 
fer the properties of the 3D fields. Their principal conclu¬ 
sions were (1) the maximum total field strength at a density 
hn > 300 cm“^ (almost entirely molecular in their sample) 
satisfies Btot.max oc hg with a = 0.65lg Qg, where we have 
estimated the errors from their Figure 4; (2) the distribution 
of field strengths is approximately uniform from very low val¬ 
ues to Btot.max; and (3) virtually all of the measured cloud 
clumps are magnetically supercritical, consistent with pre¬ 
vious results (see Section 2). Thus nature appears to favor 
the formation of magnetically supercritical clumps, which is 
a necessary condition for star formation to be possible. 
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9.1 Data Analysis 

We began by analyzing the data in a manner complementary 
to that of Crntcher et al. (2010): rather than attempting to 
infer the 3D properties of the fields, we analyzed the line-of- 
sight fields directly. As noted above, a majority of the molec¬ 
ular sources have only upper limits on their field strengths; 
by using the measured values, we are effectively stacking 
low signal-to-noise data in order to draw conclusions. This 
is valid only if the errors in the data are random and not 
systematic, but this is a prerequisite for the Bayesian anal¬ 
ysis as well. It should be noted that the fields determined 
by Zeeman observation are mass-weighted (Eq. 1), whereas 
those used in theoretical models are volume-weighted. 

1. Density dependence. Crutcher et al. (2010) found that 

the total field in molecular gas has a density-dependent 
maximum value, Btot.max = OShg®! fiG (Eq. 5). In an¬ 
alyzing their entire data set, including the HI data, they 
concluded that the distribution of held strengths was ap¬ 
proximately uniform over the range /Htot.max — Btot.max, 
with / ~ 0.03. This implies that the typical held strength 
also scales as We tested this directly by showing that 

Bios oc (uncertainty at the 95 per cent conhdence 

level), which implies that the total held scales this way as 
well. Since projection of the held onto the line of sight adds 
considerable dispersion, the dispersion in i3tot(hH) from the 
simulations is considerably less than that in i3ios(hH). In Ap¬ 
pendix E we show that the density dependence of the mag¬ 
netic held is most accurately inferred from the line-of-sight 
data by including all the data, unweighted by its accuracy. 

2. The normalized field .strength. The normalization of the 
relation between the held and the density is also signihcant. 
Since the average total held is twice the average line-of-sight 
held, (Btot) = 2(Bios) (Sec. 3), we were able to infer the 
average value of the magnetic held in the observed molecular 
clumps as a function of the density in the clumps directly 
from observation, 

(Btot) « 42 h&®| mG. (33) 

Insofar as the molecular clumps in the sample of Crutcher 
et al. (2010) are representative, this result describes the av¬ 
erage held in the molecular interstellar medium. 

3. Distribution of normalized field strengths: By compar¬ 
ing the observations of the line-of-sight helds directly with 
the projected PDFs for uniform and log-normal held dis¬ 
tributions (Appendix B) we concluded that the log-normal 
distribution is slightly better, but the difference is not sta¬ 
tistically signihcant. The primary result of Crutcher et al. 
(2010) is that the normalized total held strength Btot/h-H®® 
must have an intrinsically broad distribution; their data is 
consistent with a uniform distribution, but other distribu¬ 
tions are possible (Crutcher, private communication). The 
principal difference between the uniform and log normal dis¬ 
tributions is in the fraction of very low and very high held 
strengths: Whereas a uniform distribution has about 10 per 
cent of the clumps with helds less than 20 per cent of the 
median, a log normal has only 2.6 per cent in that range, 
and whereas a uniform distribution has no clumps with helds 
greater than twice the median, a log normal has 20 per cent 
in that range. 


9.2 Results of the Simulations 

To confront theory with observation, we have performed two 
large scale MHD simulations to study the formation of mas¬ 
sive, dense infrared dark clouds (IRDCs) (Li et al. 2015, in 
preparation). We drove Mach 10 turbulence for two-crossing 
times with no gravity; we then turned gravity on, and contin¬ 
ued driving the turbulence. For the standard linewidth-size 
relation in molecular gas, the size of the region simulated 
is 4.55 pc (Section 4). Here we have used the data from 
the gravitationally collapsing phase of these simulations to 
study the properties of the magnetic helds in the 100 most 
massive clumps identihed by CLUMPFIND in each simula¬ 
tion and have compared the results with those of Crutcher et 
al. (2010). Our two simulations have the same initial condi¬ 
tions, except that the magnetic held in the weak-held model 
is 10 times weaker than that in the strong-held model; the 
initial Alfven Mach numbers in the two simulations were 
Ma,o = 10 (weak held) and Ma,o = 1 (strong held). The 
Alfven Mach number observed in molecular clouds is typ¬ 
ically of order unity (Crutcher 1999); by the end of our 
simulations, A4a was about 0.8 and 1.6 in the strong and 
weak-held simulations, respectively. 

The differences between the observations and our simu¬ 
lations must be borne in mind. The observed data are very 
heterogeneous, coming from different molecular clouds with 
different initial conditions, size scales, and different ages, 
whereas we have carried out only two simulations that dif¬ 
fer only in the strength of the initial magnetic held. Various 
telescopes were used to obtain the observed data, whereas 
we used a single beam size of 0.025 pc to ‘observe’ our 100- 
clump sample, corresponding to the minimum beam size 
used to observe the 68 OH-I-CN sources in the sample of 
Crutcher et al. (2010). The observed clouds have masses up 
to about I 7 OOM 0 , whereas the most massive cloud in the 
simulations has a mass of about 78M©. This difference in 
masses is due both to the limited size of the simulated re¬ 
gion and to the fact that the resolution of our simulated 
beam is greater than that used for the observations of the 
most massive clouds. However, our 100-clump samples have 
clumps with densities comparable to some OH sources and 
to all the CN sources. Despite these many differences, our 
simulations are in remarkably good agreement with the ob¬ 
servations. 

1. Critical density. In Crutcher et al. (2010), the HI 
sources appear to be in a horizontal branch with Bios in¬ 
dependent of hu; the power-law relation between Btot, max 
and hn sets in above a density no — 300 cm“®, which is 
close to the critical density from a balance between turbu¬ 
lent and thermal pressure, as discussed in Section 5. The 
clumps from our two simulations also have a relation be¬ 
tween Bios and nu with horizontal branch at low density 
and a power-law branch at high density; as for the observed 
helds, the power-law branch begins at about the density at 
which the thermal pressure in the clumps balances the aver¬ 
age turbulent pressure. When gravity is hrst turned on, there 
is only a horizontal branch; the power-law branch is created 
solely as the result of gravitational collapse of clumps. 

2. Density dependence, Btot oc «)): The time-averaged 
values of a from our two simulations are 0.70 ± 0.06 and 
0.57T0.05 (uncertainties at the 95 per cent conhdence level) 
for the strong and weak-held models, respectively. By com- 
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parison, the value inferred from observation by Crutcher et 
al. (2010) is 0.65^0 05 (these uncertainties are our estimates 
based on the results presented in Crutcher et al. (2010)). The 
density scaling of the magnetic field is due to self-gravity 
since the fields were independent of density prior to turning 
gravity on. As pointed out by Crutcher et al. (2010), the 
observed value of a is very close to the value of | predicted 
by Mestel (1966) for spherical contraction of a cloud with a 
weak field. We have confirmed that the contraction of simu¬ 
lated clumps is approximately spherical. However, although 
the field in almost all the clumps is magnetically supercrit¬ 
ical, the field is not very weak: In the strong-held case it 
preserves the memory of the initial held orientation, and in 
the weak-held case it is stronger relative to gravity in the 
clumps than in the strong-held case. The B —hn relation ap¬ 
pears to be an emergent property of turbulent, magnetized, 
self-gravitating media, and a full understanding must await 
further research. We have also shown that tangling of the 
held increases with density in the weak-held case, causing 
a to be smaller for the weak-held simulation than for the 
strong-held one. As discussed in Appendix E, a systematic 
variation of noise with density can affect the inferred value 
of a. For the case at hand, in which the low-density, OH 
data is noisier than the high-density, CN data, the inferred 
value of a is smaller than the actual value by about 0.02. 
This ehect is much stronger if data less than 2 sigma is omit¬ 
ted from the sample; in that case, the inferred value of a is 
smaller than the actual value by 0.06. 

3. Mean and median values of Bios/uh®® in the simula¬ 
tions are about half and 75 per cent, respectively, of the ob¬ 
served values (Table 3). Given the many differences between 
the observed and simulated clumps, this is good agreement. 
The mean and median values remain approximately con¬ 
stant throughout the self-gravity phase of the simulations. 

4. Scaling of the magnetic field with velocity dispersion, 

a: As noted in Section 6.1, Mouschovias (1991) suggested 
that B scales with both the density and the velocity dis¬ 
persion, B (X and Basu (2000) showed that this is 

in excellent agreement with observation. We have confirmed 
that the data of Crutcher et al. (2010) are consistent with 
this scaling, although the density power law has a signif¬ 
icant uncertainty. Basu (2000) derived this relation under 
the assumption that clouds are planar and concluded that 
the agreement with observation supported this hypothesis. 
In Appendix A we show that this relation comes from an 
identity and applies to clouds of a wide range of shapes (in¬ 
cluding spherical) that are gravitationally bound and are 
close to being magnetically critical. 

5. Distribution of normalized line-of-sight field strengths. 
A K-S test shows that the median-normalized distributions 
of Bios/h-H®® of the strong- and weak-held simulation sam¬ 
ples and of the observed OH-I-CN sources are all consistent 
with being drawn from the same distribution. 

6. Distribution of normalized total field strengths in the 
simulations. The distribution of Btot/hH®®/(Btot/hH®®)nied 
from either simulation is not a uniform distribution (as con¬ 
cluded by Crutcher et al. (2010) for the observed clumps) or 
a delta function, but it is consistent with a log-normal dis¬ 
tribution. The uniform distribution of normalized total held 
strengths inferred by Crutcher et al. (2010) from observation 
implies that min(Btot, M/h-H®!) = 0.24(Bios. M/hH'®|)med- Our 
strong-held simulation has a minimum value almost 3 times 


larger. Field-line tangling increases the rms held above the 
total held; averaged over the whole clump, the minimum 
value of Brms, m/^h ,®4 is about 4 times larger than the min¬ 
imum value Btot,(i/hH®| inferred by Crutcher et al. (2010); 
the average value is about twice the minimum. 

7. The typical molecular clump in the ISM is significantly 
supercritical. The ofcsertiedmass-to-hux ratio, /i 4 .,ios, is based 
on the line-of-sight held. Bios, which is mass-weighted, and 
the column density in the central beam. For an axisymmetric 
cloud, the actual mass-to-hux ratio is based on the ratio of 
the surface density, E, in a central cylinder to the area- 
weighted held threading that cylinder in the midplane. For 
a turbulent cloud, we adopt /i<t,actual = min(E/Baiid,a'), 
where the minimum is evaluated with respect to different 
orientations and where BaUd, z' is the component of the held 
in that orientation averaged over the area of the central disc. 
We hnd 

/4<4>, actual, med — 0. 7/r4>, los, med , (34) 

substantially larger than the value (0.25 — 0.5)/i*, los, med ex¬ 
pected from projection effects, primarily due to the fact that 
the observed is based on a density-weighted held, whereas 
the actual one has no density weighting. Furthermore, the 
column densities are approximately isotropic (there is no evi¬ 
dence for signihcant battening along the held lines), possibly 
due to the turbulent nature of the held. Since the observed 
median value of /i$,ios is about 4, we infer that the actual 
median value is about 3: The typical molecular clump in the 
ISM is signihcantly supercritical. The harmonic mean of the 
actual mass-to-hux ratio is predicted to be about 1.5 based 
on our results from the strong-held model. 

8. Harmonic mean mass-to-flux ratio. The average value of 
the observed mass-to-hux ratio is based on the line-of-sight 
held, and it diverges. It is therefore better to average the 
hux-to-mass ratio, which is equivalent to using the harmonic 
mean of the mass-to-hux ratio. The same holds true for the 
actual mass-to-hux ratio if the distribution of held strengths 
extends to very low values, as suggested by Crutcher et al. 
( 2010 ). 

9. Strong vs. weak mean fields in molecular gas. The 
strong-held (initial Alfven Mach number A4a,o = 1) and 
weak-held (initial Alfven mach number Ma,o = 10) give 
substantially similar results in almost all the measurements 
discussed above because turbulence is very effective at am¬ 
plifying the initially weak held so that it is comparable to 
the held in the strong-held case, consistent with previous 
results (e.g., Lunttila et al. (2008); Padoan et al. (2010)). 
The mass-to-hux ratio based on the mean held in the tur¬ 
bulent box is 10 times higher in the weak-held case than in 
the strong-held case, but the mass-to-hux ratio in individual 
clumps is about the same in the two cases due to turbulent 
fragmentation in the weak-held case. There is a signihcant 
diherence in the accretion rates onto the clumps in the two 
simulations, with the rate in the weak-held simulation being 
slower due to the more complex held geometry. 

10. Non-ideal MHD effects. The agreement of our results 
with observations of magnetic helds in molecular cloud 
clumps is consistent with a relatively minor role for am- 
bipolar dihusion in the observed clumps. We did hnd some 
evidence for non-ideal effects in our simulations, however: 
The time-averaged value of /r# (Bj„itj_ z ') in the strong-held 
case is larger than the initial value of for the entire box. 
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This can occur only due to numerical effects. We presented 
evidence that our simulations are converged, so that our 
results are consistent with the increase being due to recon¬ 
nection in a turbulent medium (Lazarian et al. 2015). 

11. The mean field in the ISM is strong. Based on obser¬ 
vations of the alignment of polarization on large and small 
scales, Li et al. (2009) concluded that the mean field in the 
vicinity of molecular clouds is strong, corresponding to an 
Alfven Mach number Ma{{B)) < 1. Our simulations sup¬ 
port this conclusion: There is negligible correlation between 
the orientation of the clump fields with the mean field in the 
weak-held simulation, whereas the correlation in the strong- 
held simulation is comparable to the observed values. 
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APPENDIX A: THEORY: THE MAGNETIC 
FIELD VS. GRAVITY 

The relative importance of magnetic fields and gravity is 
usually parametrized by the ratio of the mass to the mag¬ 
netic critical mass (Mouschovias & Spitzer 1976), 

Mfl, = c<s. . (Al) 

where <!> is the magnetic flux and c<s> is a numerical con¬ 
stant. The value of c$ depends on whether M<j is for the 
entire cloud or for the central flux tube; in the former case, 
it also depends on the flux distribution. The flux distribu¬ 
tion can be described in terms of the equilibrium shape of 
the cloud in the absence of gravity, so that the field is uni¬ 
form. Mouschovias & Spitzer (1976) considered the flux dis¬ 
tribution corresponding to a cloud that is spherical in the 


absence of gravity; Tomisaka et al. (1988) considered sev¬ 
eral cases, including both a spherical and a cylindrical ini¬ 
tial cloud. For the spherical case and taking 4> to be the flux 
threading the entire cloud, Mouschovias & Spitzer (1976) 
found c<s> = 0.126. Subsequent calculations by Tomisaka 
et al. (1988) gave ~ 0.12 for this case. On the other 
hand, if one focuses on the central flux tube, Tomisaka et 
al. (1988) showed that c<j = 0.17 — 0.18 for the three cases 
they considered, so that it is relatively independent of the 
flux distribution. This value is very close to the value ap¬ 
propriate for a slab threaded by a uniform, perpendicular 
held, c<j = l/27r (Nakano & Nakamura 1978). Since this is 
very close to the value for the central hux tube and since 
Crutcher et al. (2010) used this value, we shall adopt it also, 
both for the clumps that form in our simulation and for 
the simulation box (Section 4) as a whole. The normalized 
mass-to-hux ratio (henceforth, we shall refer to this as just 
the ‘mass-to-hux ratio’) is then 


_ Af ^1/2 S 


3.80 


Ah, 21 


(A2) 


where S is the mass surface density. Ah, 21 is the column 
density of H nuclei in units of 10^^ cm“^, B is the magnitude 
of the magnetic held, and — B/{I /rG). 

It is convenient to relate /i* to other dimensionless pa¬ 
rameters that characterize the system. Since we wish to do 
this for the simulation box as well as in clumps, both simu¬ 
lated and observed, we hrst dehne two geometrical factors. 
For the simulation box, let i = io he the box size, whereas 
for a spherical cloud, let i be the diameter. We then dehne 


A = CAf, (A3) 

V = cvC, (A4) 


where ca = cv ~ 1 for a box; for a sphere, ca = 7r/4 
and Cv = rr/6; and for an ellipsoid that is axisymmetric 
about the 2: axis with semi-major axes R = 112 and Z, 
Ca =1^1^ for fhs surface normal to 2 and cv = (ti/ e))Z/R. 
For example, the surface density, S, of a cloud is related to 
its mean density, p, by E = M/A = (cy/CA)pt. 

Next, we compare the kinetic and gravitational energies 
of a cloud. Let a^t be the mass-weighted, non-thermal ID ve¬ 
locity dispersion, Vrms = \/3crnt the 3D non-thermal velocity 
dispersion, Cs the isothermal sound speed, and M = Urms/cs 
the 3D Mach number. The mass-weighted, total ID velocity 
dispersion, am, is then given by 

<^1D =C^ + CTnt = (f + = CMalf (A5) 

Since the motions in molecular clouds are often highly su¬ 
personic, CM is generally of order unity. To measure the 
relative importance of kinetic energy (including thermal mo¬ 
tions) and self-gravity, we use the virial parameter (Bertoldi 
& McKee 1992) 


_ 5amt _ PcMantt 
~ 2GM ^ 2GM 


(A6) 


where we have replaced the radius in the original definition 
by 1/2. For centrally concentrated clouds, Ovir is somewhat 
greater than twice the ratio of the kinetic energy to the 
gravitational energy (Bertoldi & McKee 1992). The mass- 
to-flux ratio can then be expressed as 
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57VCvCm 

6c^ftvir 


1/2 

Ma, 


(A7) 


where A4a = Wrms/wA is the 3D Alfven Mach number and 
va ~ B/ {4:TTpy^^ is the Alfven velocity. For a sphere, the 
nnmerical factor is (STrcv/hc^)^^^ = (20/9)^^^ = 1.49, only 
slightly less than the value for a box, = 1.62. 

Furthermore, both these factors are of order unity: It follows 
that highly supersonic clouds with oivir ~ Ma ~ 1 are close 
to magnetically critical, /i^. ~ 1. In Section 7 we evaluated 
with several different values of B, and it must be kept 
in mind that the same value of B must be used in both 
and A4a when applying Equation (A7). 

Note that solving Equation (A7) for B gives a relation 
of the form derived by Mouschovias (1991) and Basu (2000) 
for planar clouds. 


/ IOtt^Cv 

y C^Qvir 


1/2 -1/2 
P UiD 

P<S> 


(AS) 


Note that this equation is actually an identity; it follows 
directly from the dehnitions of the quantities entering it. 
Physics enters when it is assumed that the cloud is self- 
gravitating ((Tvir — 1) and magnetically critical (/i$ ~ 1), as 
Mouschovias (1991) did. Basu (2000) argued that the fact 
that this relation provided a good ht to the data suggested 
that clouds are indeed flattened; however, we see that a re¬ 
lation of the same form applies to spherical clouds as well. 


A1 Relation between BjriH^ and the Magnetic 
Critical Mass Mb 


An alternative form for the magnetic critical mass is 

M M<s> 

Mb = ^ = 


P% 


9 ’ 

ni 


(A9) 


which was found by Mestel & Spitzer (1956) and related to 
Mii by Mouschovias & Spitzer (1976). Note that since M oc 
P and pij, oc M/<1> oc £, the critical mass Mb is independent 
oi£: 

/ 3 \ f>3 

Mb = 2.15 X 10-^ ^ Mq. (AlO) 

\Cv J » 1 ' h .4 


For an ellipsoid we have c\lcy = (97r/16)(i?/Z)^. Equations 
(A9) and (AlO) imply 


- 2/3 


16.7 


I caP<s> 



(All) 


Generalizing the scaling relations (15) and (16) given in Sec¬ 
tion 4.1 to include thermal motions (cm > 1), this becomes 


B, 

- 2/3 


= 11.3 




(A12) 


Thus, the ratio of the magnetic held to a power of the den¬ 
sity, which is dimensional, is determined by the temperature, 
the coefficient in the linewidth-size relation and the dimen¬ 
sionless parameters that characterize the problem. 

It is important to note that the observed mass-to-flux 
ratio is p<s>{Bios), where Bios is the mass-weighted held along 
the line of sight, whereas the value used in theoretical calcu¬ 
lations, which assume axisymmetry, is p<s>{Bp), where Bp is 


the poloidal held threading the midplane of the entire cloud 
(Mouschovias & Spitzer 1976) or the central region of the 
cloud (Tomisaka et al. 1988) and is not density-weighted. 


APPENDIX B: PROJECTED FIELD 
DISTRIBUTION 


The distribution of the observed values of Bios,is 
shown in Fig. 2. In order to relate a given model for the total 
helds to the observed line-of-sight helds, we must determine 
the probability distribution of the line-of-sight held, ip{Bios), 
from that of the total held, <j){Btot)- This has been done by 
Heiles & Crutcher (2005), 

^{Bios)= f (Bl) 

JbIos 

One can readily show that this implies that 

(Aos) = -^{B^oth (B2) 

n -I- 1 

in particular, (Bios) = |(Btot). Equation (B5) can be diher- 
entiated to give 


4'{Btot) 


d'i)i(Bios) 
d In Bios 


^los — -®tot 


(B3) 


but this is not very useful in analyzing data because the 
derivative amplihes the noise. In terms of the normalized 
held, b, dehned in Equation (7), Equation (B5) becomes 


'Ipibloa) = 



dbtot . 


(B4) 


In order to relate a given model for the total helds to the 
observed line-of-sight helds, we must determine the probabil¬ 
ity distribution of the line-of-sight held, ip{Bios), from that 
of the total held, (f>{Btot)- This has been done by Heiles & 
Crutcher (2005), 


^(Blos) 


r 4>{Btot) 

JbIos 


dBtot. 


(B5) 


One can readily show that this implies that 

(Aos) = ^(Brot); (B6) 

n -I- 1 


in particular, (Bios) = 5 (Btot). Equation (B5) can be diher- 
entiated to give 


4’{Btot) 


dip{Bios) 
d In Bios 


^los — -®tot 


(B7) 


but this is not very useful in analyzing data because the 
derivative amplihes the noise. In terms of the normalized 
held, b, dehned in Equation (7), Equation (B5) becomes 


'Ipibloa) = 



dbtot . 


(B8) 


We shall consider three possible distributions of held 
strengths. The hrst is the uniform distribution of Crutcher 
et al. (2010), 

<j>{btot) J- TTT 

J- J )0tot, max 

[H(b tot y^tot,max) H(6tot ^tot,max )],(B9) 
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Table Bl. Models of the Observed Fields 



Observation 


Model 




Uniform 

Log normal 

Exponential 



/ = 0.007 

CT = 0.825 


{Bios,m/«&,®|) 

22.45 

19.91 

24.68 

22.69 

4 )med = ^los, med 

11.94 

14.93 

15.01 

15.73 

(Bios, fj./^\os, med 

1.88 

1.33 

1.64 

1.44 

= 2{Sios,mK-6|> 

44.90 

39.81 

49.37 

45.38 

(Btot,M/^^&,®|)med 

- 

39.81 

35.15 

38.10 

(Btot, /i/'^^'^ 4 )med/^los, med 

- 

2.67 

2.34 

2.42 

fraction of (Btot, < ^los,med 

- 

18% 

15% 

15% 

fraction of (Btot, > 3foios,med 

- 

44% 

38% 

39% 

Reduced of fit 

- 

1.85 

1.45 

1.35 


where H(a;) is the Heaviside step function. For such a dis¬ 
tribution, the distribution of the line-of-sight field (Eq. B5) 


t/’(6los) = 


1 + / 


(1 - /)4(6io.) 


In min 


4(&los) 1 

(l + /)6io=’ / 


(BIO) 


for feios ^ &tot,max and t/;(6ios) = 0 otherwise. In expressing 
this result, we have taken advantage of the general result 
(Bios) = |(Btot), which implies 

u _ 4(&los) ,., 3 ^ ., ^ 

Otot,max — ^ _j_ y ' 

The mean and median of the total density-normalized held 
are 

(^tot) — (^tot)med ~ /)^tot, max. (B12) 

By contrast, the median of the line-of-sight held is less than 
the mean; for / = 0, it is (&ios)med = 0.748(6ios)- Since 
(btot) = 2(6ios), and since (btot) = (6tot)med, it follows that 
for / = 0 the median line-of-sight held is related to the 
median total held by (&ios)med = 0.374(6tot)med- 

A second possible form for the distribution of total held 
strengths is a log normal, for which the distribution of held 
strengths is given by 


0(l^tot)d6tot — 


(27r)i/2, 


■ exp 


{x+^a'^y 

2a2 


dx. 


(B13) 


a:ios = In 


^los 

(^tot) 


T ^lOS 

in —t;— r- 

2(&los) 


(B17) 


The median value of &tot is related to the average value by 
&tot,med = (6tot) exp(—cr^/2). There does not appear to be 
a simple expression for the median value of the line-of-sight 
held. However, we can obtain an approximation as follows. 
For (T —>■ 0, the distribution of btot is a delta function and 
-!/>(&ios) is uniform between 0 and (btot) = btot’, as a result, 
(&ios)med = |6tot- On the other hand, in the limit cr ^ 1, the 
dispersion in btot dominates the dispersion due to the projec¬ 
tion onto the line of sight, and (feios)med = (f>ios)e~'^ just 
as in the case of btot, meu and (btot) (this can also be shown 
by direct calculation). An approximation that reduces to 
these two limits is (6ios)med — (6ios)e“‘^ It follows that 
the dispersion of the total held can be inferred directly from 
observations of the line-of-sight held from 


(T ~ 2 In 


{bios') 


(bit 


(B18) 


The third distribution we shall consider is an exponen¬ 
tial distribution of the line-of-sight held. 

Equation (B7) then implies that this corresponds to a dis¬ 
tribution of the total held. 


where 


* = ln^. (B14) 

\0tot) 

Note that 0 is not a log normal itself, but rather a log normal 
times dx/dbtot ~ l/&tot. The distribution of the line-of-sight 
held (Eq. B5) is then 


'ipibios) 


(8'7r)i/2(j, 

(^los) 

poo 

{x+icy 

/ exp 

2cr2 


dx, 


4(&los) 


erfc 


Xlo, 


+ IC 


^2(7 


(B15) 

(B16) 


where erfc(a:) is the complementary error function and 


This distribution differs from the hrst two in that it is not 
symmetric. 

Model hts to the observed line-of-sight data from 
Crutcher et al. (2010) are summarized in Table Bl. We eval¬ 
uated the goodness of ht with a test in log(i3ios, 

For the uniform distribution, = 1.85. The best-ht value 
oi f — min(Btot/Btot, max) is 0.007, smaller than the value 
of 0.03 found by Crutcher et al. (2010), which included HI 
sources as well. The reduced x^ values of the lognormal and 
exponential distribution (at the bottom of the table) are 
slightly smaller than that of the uniform distribution, im¬ 
plying that they provide slightly better fits. Table Bl also 
gives several characteristic values of Bios for the three theo¬ 
retical models, as well as for the observations summarized in 
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Crutcher et al. (2010). Both the uniform and lognormal dis¬ 
tributions have means that differ by about 10 per cent from 
the observed mean value, whereas the exponential function 
has a mean that is significantly closer. The median values of 
all three models are larger than the observed median value 
by about 25 per cent. The ratio of the mean to median, 
which is a measure of the shape of the distribution, for the 
lognormal distribution is closer to the observed value than 
for the other two distributions. 


APPENDIX C: SELECTION REQUIREMENTS 
OF THE 100-CLUMP SAMPLE 

To select clumps from simulations for comparison with ob¬ 
servations, we impose several requirements: 

1. Resolution: The 68 OH and CN Zeeman observations in 
Crutcher et al. (2010) employed several different radio tele¬ 
scopes. There is information on the beam sizes for the 34 OH 
Zeeman observations using the Arecibo telescope (Troland 
& Crutcher 2008) and the 14 CN Zeeman observations us¬ 
ing the IRAM-30 m telescope (Falgarone et al. 2008). It is 
important to note that Bios is measured in a single telescope 
beam (see Section 2), and for resolved clumps this is smaller 
than the whole clump, which is usually defined as the re¬ 
gion inside the half-maximum contour in the intensity map. 
Higher resolution observations will thus probe the dense in¬ 
ner region of a clump better than observations with a larger 
beam. The observed clumps we compare with have a mini¬ 
mum radius of 0.03 pc. The minimum beam size of the 48 
OH and CN clumps is 0.025pc. The remaining OH sources 
are most likely observed with a larger beam size. Taking into 
account the projection effect from observations, the 3D dis¬ 
tance between two clumps of minimum projected distance 
will be « 0.058 pc. We therefore merge all clumps that are 
separated by 0.06 pc or less, since they would be seen as 
single clumps in the observations. A possible correlation be¬ 
tween the clump properties and the beam size is discussed 
in Section 6.1. 

2. Density: In Fig. 1 of Crutcher et al. (2010), the ob¬ 
served OH and CN clumps appear to have a power-law re¬ 
lation between the maximum Bios and the density, whereas 
the HI clumps are on a horizontal branch, with the max¬ 
imum Bios independent of density. The power-law relation 
starts roughly at hn ~ 300 cm“^. As discussed in Section 
2, this density corresponds approximately to the critical 
density, nn.crit, at which the thermal pressure of the gas 
balances the turbulent pressure of the ambient interstellar 
medium. Our models have a much higher turbulent pressure, 
Tturb/fc = 1.92 X 10® K cm“®, because we do not have suf¬ 
ficient dynamic range to include the diffuse ISM as well as 
dense CN clumps. At T = 10 K and Xtot = 0.6 for molecular 
gas, this implies nu, crit = 3.2 x 10"^ cm“®. Since the focus 
of our research is on magnetic fields in molecular gas, all of 
which is above the critical density, we therefore require our 
clumps to have minimum density larger than 4 x 10'* cm“®. 
Below, we shall see that this density corresponds approx¬ 
imately to the beginning of a power-law relation between 
Bios and density, just as in the observations. Furthermore, 
these high densities are a result of self-gravity; only a few 
clumps have densities above this threshold prior to the time 
at which self-gravity is turned on. 


3. Mass: There are still a few hundred clumps after filter¬ 
ing out the clumps that are too close to adjacent clumps 
or are too low in mean density. On the other hand, the ob¬ 
served sample of molecular clumps with fields measured by 
the Zeeman effect is less than 100. In order to have a sam¬ 
ple that is comparable to the observed one, but has enough 
clumps to obtain reasonably accurate measured mean val¬ 
ues, we choose the 100 most massive of the simulated clumps 
that survive the filtering process. All these clumps have at 
least 10000 grid cells, so they are well-resolved; fortuitously, 
they also have sizes B > 0.03 pc, the minimum resolution of 
the observations that we are comparing with. Note that our 
clumps are less massive than the observed ones: The total 
mass in the simulated region is only 3110Mq and not many 
of the simulated clumps have masses larger than IOMq, 
whereas the most massive clumps in the observed sample 
exceed 10®M©. We also looked at the 100 clumps with the 
largest volume, 40 per cent of which were not in the most 
massive sample, and found that (&tot), {bios), and (/i*(i3ios)) 
in the two samples agree to within approximately 1 sigma. 

In order to study the properties of the clumps in the 
absence of self-gravity, we also consider non-self-gravitating 
100-clump samples from the strong and weak field simula¬ 
tions, in which there is no density threshold. The clumps in 
these samples have densities hn ^ 2 x 10'* cm“®, about half 
the minimum in the self-gravitating case. 


APPENDIX D: 100-CLUMP SAMPLE 
CONVERGENCE STUDY 

We compare magnetic field properties of clumps obtained 
from CLUMPFIND in this paper with observations. There¬ 
fore, our convergence study focuses on the magnetic field 
properties of clumps for the whole clump, the central region 
of the clumps, and the central region of the clump convolved 
with a Gaussian beam. We compare the clumps from the 2- 
level refinement, strong-held model with a 3-level refinement 
run with exactly the same initial conditions. With our com¬ 
puting resources, we can afford to run the 3-level model only 
up to 0.4tff, when the power-law branch of the clumps is well 
developed; we therefore compare the two runs at that time. 
See Section 5 on how the 100-clump samples are obtained 
using CLUMPFIND and Section 6 on how the size of the 
convolved Gaussian beam region is chosen. 

Table A.l shows the convergence results for the 100- 
clump sample convolved with a Gaussian beam of 0.025 pc 
radius. Only line-of-sight quantities and quantities, such as 
Btot that are computed based on line-of-sight quantities, are 
included in this table. The uncertainties given in each case 
are the standard error of the mean; all quantities are con¬ 
verged to within this value. Table A.2 shows the means and 
medians of the volume- and mass-weighted magnetic field 
strength, the mass-to-flux ratio, the column density, and the 
Alfvenic Mach number, both for whole clumps and for small 
central cylinders of radius 0.025 pc. We use a cylinder rather 
than a Gaussian beam here since we are testing for conver¬ 
gence of the actual mass-to-flux ratio, not the observed one. 
Almost all quantities from the two runs are converged within 
1 sigma, except the mass-weighted Srms.M is slightly larger 
than 1 sigma but well within the sum of the two errors of 
the mean. With the same sample size for both the 2-level 
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Table A.l. Convergence Study of 100-Clump Sample Convolved 
With a 0.025 pc Gaussian Beam at 0.4 tfj 


Model" 

2-level 

3-level 

(Bios)" 

76.0 ± 7.4 

70.4 ± 7.4 

(-^los)med 

43.6 

42.4 

(Stot> 

151.6 ± 21.1 

148.2 ±20.6 

(-Stot)med 

95.7 

90.9 

(Bios/4®^)med 

0.028 ± 0.003 

0.028 ± 0.003 

0.023 

0.021 


0.054 ± 0.003 

0.055 ± 0.004 

(Btot ) rned 

0.049 

0.045 

(M<i>,los (Bios)} h 

2.0 ±0.1 

2.0 ±0.2 

(M<i>,los (-®los))med 

2.6 

2.7 

(Nh)^ 

4.04 ± 0.29 

4.05 ±0.23 

(A^H)med^ 

2.61 

2.74 

(^H,los) 

1.91 ± 0.13 

1.89 ±0.13 

(^H,los)med^ 

1.10 

1.11 


“ Two models labeled by the number of refinement levels. 

All magnetic field strengths are in units of fiG. 

Harmonic mean of the line-of-sight mass-to-fiux ratio is 
computed using Bios- 

Column density in unit of 10^^ cm“^. 

® Volumn density computed from Nii/{2 X rdump) in unit of 10® 
cm“®. 


and 3-level 100-Clump samples, the difference is not statis¬ 
tically significant. Therefore, the 2-level refinement run is 
converged at 0.025 pc radius size scale. 


APPENDIX E: SAMPLING THE 
DISTRIBUTION OF THE LINE-OF-SIGHT 
MAGNETIC FIELD 

The OH and CN sources discussed in this paper are 
from three separate papers by Crutcher (1999), Troland 
& Crutcher (2008), and Falgarone et al. (2008), using dif¬ 
ferent telescopes at different wavelengths. From Fig. 1 in 
Crutcher et al. (2010), the measured line-of-sight magnetic 
field strength Bios ranges from a few to a few thousand /rC. 
Forty of the 68 OH-I-CN sources have measured values of 
Bios < 2ct uncertainty. Because the noise level is so large, 
the true values of Bios could be significantly different than 
the values listed in Table 1 in Crutcher et al. (2010), which 
could have a significant effect on the inferred the power-law 
index a that relates Btot and nn- Here we use a numerical 
experiment to determine how to infer the value of a from 
noisy data. 

In this experiment, we assume that the normalized 
total magnetic fields, Btot/h-H, have a uniform distribu¬ 
tion between 0 and Btot, max/uH^®, close to the distribu¬ 
tion inferred by Crutcher et al. (2010). We then generate 
an 8000-point sample of from this distribution. 

Since the observed OH sources have densities in the range 
nn = 1 X 10® — 2 X 10"^ cm“®, and the CN sources are in the 
range nn = 2 x 10® — 4 x 10® cm“®, these 8000 points are 
generated in two equal groups in these two density ranges. 
We then randomly project Btot to the line-of-sight to ob¬ 
tain 8000 points in the Bios — hu plane. A least squares fit 
of log(i3ios) vs. log(nH) gives a — 0.65±0.01, where it should 


Table A. 2. Convergence Study of 100-Clump Sample of the 
Whole Clump and Inside the Central 0.025 pc region at 0.4 tg 


Model" 


2 -level 

3-level 

Clump Center 

(Blos)^ 

70.4 ± 6.2 

67.7 ± 7.0 


(-^los)med 

41.8 

43.3 


(-^los)) 

2.1 ± 0.1 

2.0 ±0.1 


(M^,1os (-^los))med 

2.8 

2.7 


(A^h>^ 

4.64 ± 0.29 

4.46 ±0.29 


(-^H)med 

2.84 

2.78 


(l(B>vl> 

83.7 ± 3.1 

83.5 ±3.1 


(l(B>vl)med 

79.2 

77.7 


(Brms,V ) 

117.6 ±4.5 

121.3 ±4.8 


(BymSfV )med 

113.3 

114.3 


(-Btot) 

143.0 ± 17.0 

141.3 ± 19.6 


(Btot)med 

93.7 

92.7 


(Brms,M}^ 

223.5 ±31.5 

270.4 ± 41.5 


(Bxms,M )med 

129.9 

127.6 

Whole Clump 

(l(B>y|) 

68.0 ± 2.4 

67.7 ± 2.3 


(l(B>y|)nsed 

65.0 

66.6 


(Brms, V ) 

79.8 ± 2.5 

78.4 ± 2.4 


(Bxxas^V )med 

78.2 

77.5 


(-Btot) 

108.7 ± 11.7 

98.6 ±8.6 


(Btot)med 

78.5 

76.2 


(Brms,M ) 

154.4 ± 20.7 

165.2 ±23.8 


(Brms,M )med 

96.5 

91.9 


(/i$(Btot, v))h^ 

1.70 ± 0.09 

1.62 ±0.09 


(/i<|) (B^ot, V ))med 

1.83 

1.77 


(/i<|> (|B|tot, V ))med 

1.75 

1.68 


(Ma) 

0.99 ±0.05 

1.01 ± 0.05 


(A^A)med 

0.86 

0.86 


“ The two models are labeled by the number of refinement levels. 
® All magnetic field strengths are in units of pG. 

The line-of-sight mass-to-fiux ratio is computed using Bios. 
Column density in units of 10®® cm“®. 

® We have added the subscript ‘M’ to the mass-weighted rms 
field to emphasize the distinction from Bj-^os V i which appears in 
the main text. 

■f The harmonic mean mass-to-fiux ratio for whole clumps. For 
each clump, the mass-to-fiux ratio is determined along 1152 
different orientations using the volume weighted field, Btot Vt 
and the minimum value is adopted. (Section 7.2) 


be recalled that 0.01 is 2a. This confirms that we have cor¬ 
rectly implemented the model of Crutcher et al. (2010). 

Next, we select 40 random samples of Bios, each with 
the same number of low-density points (54) and high-density 
points (14) as in the observed data. These points represent 
the ‘true’ values of this quantity. We then add noise to these 
samples. Since the CN data has less noise than the OH 
data—8 of the 14 CN sources are greater than 2a, whereas 
only 20 of the 54 OH sources are—we determine the values 
of the standard deviation a for the high- and low-density 
groups separately so that each sample has the same number 
of values of Bios above 2a as in the data. We then determine 
the ‘observed’ Bios for each data point by randomly draw¬ 
ing from a Gaussian distribution, using the ‘true’ Bios and 
a as the mean and standard deviation, respectively. With 
the noise, the ‘observed’ value of Bios can be negative; when 
this happens we retain only the magnitude of Bios since we 
are interested in the distribution of field magnitudes. 
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Once we have the 68 data points including noise for 
each of the 40 samples, we determine a in each sample from 
a least squares fit of log(-Bios) vs. log(hH) in which we treat 
all the data points equally. The result is a = 0.63 ±0.13; the 
standard error of the mean is 0.01. This value is close to, 
but slightly less than, the value a — 0.65 used in generating 
the initial 8000-point data set. The decrease in a is to be 
expected: since the low density points are noisier than the 
high density ones, the negative values are larger, as are the 
corresponding magnitudes. As a result, noise increases the 
average magnitude of the low density points more than that 
of the high density ones and the slope is reduced. Although 
we have been careful to reproduce the sampling of the data, 
that does not affect the value of a: Adding noise to the 8000 
points in the same manner as we did for the samples (i.e., so 
that 34/54 of the 4000 low density points are less than 2a 
and 6/14 of the 4000 high density ones are), a least squares 
fit gives a — 0.63 ± 0.01, just as before. We conclude that 
the greater noise in the low-density data reduces the inferred 
value of a by about 0.02 below the true value. 

The effect of the different accuracies of the low and high 
density data can also be seen if we exclude data points that 
are < 2a; we then find a is 0.59 ± 0.05, more than 2a less 
than the correct value of 0.65. Elimination of the data below 
2a increases the average value of the retained OH data more 
than of the retained CN data, thereby reducing the slope. 
Furthermore, this value of a is more than la less than the 
value obtained by treating all points equally. We conclude 
that it is more accurate to treat all the points equally in the 
least squares fitting than to omit the points that are less 
than 2a. 
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